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^~H Abstract. We consider the stochastic evolution of a (1 + 1) -dimensional polymer in the de- 

^^ pinned regime. At equilibrium the system exhibits a double well structure: the polymer lies 

vN (essentially) either above or below the repulsive line. As a consequence one expects a metastable 

r — H behavior with rare jumps between the two phases combined with a fast thermalization inside 

I ^ each phase. However the energy barrier between these two phases is only logarithmic in the sys- 

tem size L and therefore the two relevant time scales are only polynomial in L with no clear-cut 
separation between them. The whole evolution is governed by a subtle competition between the 
diffusive behavior inside one phase and the jumps across the energy barriers. In particular the 
^ usual scenario in which the tunneling time coincides with the exponential of the energy barrier 

^^ breaks down. Our main results are: (i) a proof that the mixing time of the system lies between 

Q, L2 and L2+^; (ii) the identification of two regions associated with the positive and negative 

phase of the polymer together with the proof of the asymptotic exponentiality of the tunneling 
time between them with rate equal to a half of the spectral gap. 
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o 

■^ 1. Introduction, model and results 

^^ Random polymers are commonly used in statistical mechanics to model a variety of inter- 

<^ esting physical phenomena. A rich class of models with a non trivial behavior is obtained by 

^""1 considering a simple random walk path interacting with a defect line in the thermodynamic 

J> limit when the length of the path tends to infinity. The equilibrium of these so-called polymer 

pinning models has been studied in depth in the mathematical literature, and the associated 

localization/delocalization phase transition is, nowadays, a well understood phenomenon, even 

C^ in the presence of non homogeneous interactions; see [9] for a recent survey. 

Markovian stochastic dynamics of random pinned polymers, on the other hand, have received 
much less attention from a mathematical point of view. Besides their importance in bio-physical 
applications (see e.g. [5, 6] and references therein), the stochastic evolution of polymer models 
poses new challenging probabilistic problems from many points of view and the connection 
between the equilibrium and dynamical properties of the model is still largely unexplored. In 
particular, we feel that the problem of how the polymer relaxes to the stationary distribution 
(time scales, overcoming of energy barriers, metastability, patterns leading to equilibrium) still 
lacks a satisfactory solution even in the simplest homogeneous models; see [7] for some initial 
results in this direction. 

In this paper we consider the dynamics of a homogeneous polymer model interacting with a 
repulsive defect line with two main motivations in mind: 
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(i) the repulsive regime is characterized by a relaxation to equilibrium occurring on a time 
scale certainly much larger [7] than the usual diffusive one which is typical of the e.g. 
neutral case^ [21]. The new scale is clearly the result of a complicate competition in the 
polymer evolution between diffusive behavior and jumps against energy barriers, 
(ii) The whole relaxation mechanism should show certain typical features of metastable evolu- 
tion but in a very a-typical context^ in which the relevant relaxation time scales are only 
polynomial in the size L of the system (i.e. the energy barriers are only logarithmic in 
L), with little separation between the mixing time inside one phase and the global mixing 
time . A signature of this fact can be found in the anomalous growth with L of the global 
mixing time, a growth which is much more rapid than the naive guess based on the usual 
rule Tinix ~ exp(A£'), with AE the so called activation energy. In order to appreciate the 
novelty of such a situation it is useful to compare it to another well known case, namely the 
Glauber dynamics for the low temperature Ising model in a square box with free bound- 
ary [14], for which a very precise analysis of the metastable behavior was possible exactly 
because of a sharp separation, at an exponential level, between the two time scales. 

1.1. Dynamics of the polymer pinning model. Let Q = 0,2l denote the set of all lattice 
paths (polymers) starting at and ending at after 2L steps, L G N: 

n = {r^e 1?^+^ : r]_L = r]L = 0, r?^+i = r/^ ± 1 , x = -L, . . . ,L - 1} . 

The stochastic dynamics is defined by the natural spin-flip continuous time Markov chain with 
state space Q. Namely, sites x = — L + 1,...,L — 1 are equipped with independent rate 1 
Poisson clocks. When site x rings, the height r]x of the polymer at x is updated according to 
the rules: if rjx-i = 'Hx+i i 2 then do nothing; if r]x-i = "Hx+i = h, and \h\ ^ 1, then set 
Tjx = h ziz 1 with equal probabilities; if rjx-i = Vx+i = +1, then set rjx = with probability 
^nr[ and rjx = 2 with probability jq-j-; similarly, if rjx-i = t]x+i = — 1, then set rjx = with 
probability -^—^ and r]x = —2 with probability in:^. Here A > is a parameter describing the 
strength of the attraction (A > 1) or repulsion (A < 1) between the polymer and the line h = 0. 
The infinitesimal generator of the Markov chain is given by 

L-l L-l 

mv)= E ^-,+(^)[/(^''^)-/(^)]+ E r^Av)[fiv''n-m], (i-i) 

x=-L+l x=-L+l 

where: / is a function ^2 i— )■ M; rj^' denotes the configuration which coincides with r] at every 
site y / X and equals rjx ± 2 at site x; the rates rx^± are zero unless tJx-i = Vx+i, and if 
Tjx-i = rix+i = h they satisfy rx^± = \ ioi h ^ ±1, and Vx,^ = x+T = 1 " 'rx,±j for h = ±1. 

The process defined above is the heat hath dynamics for the homogeneous polymer pinning 
model, with equilibrium measure vr = it2l on Q. defined by 

4l(^) = ^^, (1-2) 

^2L 



^In the neutral case (absence of an interaction between the polymer and the line) our process is nothing but 
the usual (finite) symmetric simple exclusion model. 

^Over the years there have many different formulations of "metastability" . We refer in particular to [18, 20, 4] 
and to the recent contributions [2, 3] where, as in our case, energy barriers are only logarithmic in the characteristic 
size of the system. We feel however that our situation does not fit completely in any of the quoted contexts. 

We cannot necessarily be very precise here about the exact meaning of these mixing times but their definition 
and their role will appear clearly later on. 
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where N{ri) = #{x G {—L + 1, . . . , L — 1} : r]x = 0} denotes the number of zeros in the path 
r/ G $7 and Z^^ = J^rt'en"^^ • ^'^^ every A > 0, L G N, vr = 7r2^ is the unique reversible 
invariant measure for the Markov chain. 

1.2. Relaxation to equilibrium. The equihbrium properties of the polymer pinning model 
have been studied in detail, cf. e.g. [10] or [9, Section 2] for an extensive review. In particular, 
it is well known that, under the measure 7r22,, for A > 1 the path is strongly localized with a non 
vanishing density of zeros, while for A < 1 the path is delocalized with vL height fluctuations 
and with the number of zeros stochastically dominated by a geometric random variable with A- 
dependent parameter. The dynamical counterpart of this localization/delocalization transition 
has not been fully understood yet. Some progress in this direction has been reported in [7], 
where various bounds on the spectral gap and mixing time of the Markov chain were obtained, 
together with estimates on the decay of time correlations. 

We recall that the spectral gap is the smallest nonzero eigenvalue of —C, and one is often 
interested in the relaxation time Tj-ei = 1/gap which governs decay to equilibrium in L^(7r), while 
the mixing time T^[x{d), for S G (0, 1), is the smallest time t such that 

max ||P((ry, •) — 7r|| ^ (5 , (1-3) 

where Pt{ri, ■) denotes the distribution of the Markov chain at time t with initial state r] £ Q, and 
11^ — 1^1 1 is the usual total variation distance between two probability measures. When <J = ^ we 
often write simply Tmix instead of Tinix((5). With these conventions one has Tj-ei ^ Tmix always, 
and the inequality is strict in general. 

A dynamical phase transition occurs when we move from the localized regime A > 1 to the 
delocalized regime A < 1. It was shown in [7], see Theorem 3.4 and Theorem 3.5 there, that for 
A ^ 1 one has Tj-d = 0{L^) and Tmix = 0{L^ log-L) for all L, while for A < 1 one has 

r,ei^Li-S (1.4) 

for all e > 0, provided L is large enough. 

We refer to [7] for results and conjectures concerning the localized regime A > 1. Here, we 
consider the delocalized regime, i.e. in the repulsive case A < 1. The first question we address 
concerns an upper bound on the relaxation time Tj-d and the mixing time Tmix- It is worth 
noting that even a crude polynomial bound is non-trivial. We refer to [15, 7] for polynomial 
bounds for the model with a horizontal wall at level zero, i.e. when lattice paths are constrained 
to be non negative. On the other hand, without the wall constraint, the equilibrium measure vr 
is known to be concentrated, as L — t- cx), on configurations in which the density of monomers 
in the upper (lower) half plane is approximately either one or zero. However a mathematical 
working definition of a candidate for the plus or minus phase for the polymer is not so obvious 
and we have been forced to introduce a mesoscopic parameter £ {i.e. L ^ £ ^ 1) and define Q,^ 
by 

n+ = {r]£n : T]x>o, -L + i<x<L-i}, n~ = -n+ , (1.5) 

where, for any set A of polymer configurations —A = {rj G Q : —rj G ^}). 

The presence of the two phases associated to $7 dramatically changes the relaxation scenario, 
with a bottleneck at the set O \ (fi^ U Q^). As explained in [7, Section 6], one may suspect that 
Trei ~ La is the correct asymptotic behavior in the delocalized regime. Let us briefly recall the 
heuristic reasoning behind this prediction. 

The time to reach equilibrium can be roughly thought of as the time needed to switch from, 
say, 0~ to Q'^ . A point x such that t]x = and rjx^i ^ "Hx+i is called a crossing of the polymer. 
Note that any zero (and therefore any crossing) x must belong to the set El of points in the 
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segment {—L + 2, . . . , L — 2} which have the same parity as L. Since there are typicahy very few 
zeros at equiUbrium, one may consider the extreme case where at most one crossing ^ is ahowed 
at ah times. In this case, the time evolution of ^ should be essentially described by a suitable 
birth and death process or random walk on El; see Figure 1. 






Figure 1. From left to right a snapshot sequence of the motion of a single 
crossing S, which allows the system to switch from a mostly negative to a mostly 
positive configuration. 

From equilibrium considerations, one knows that this random walk should have reversible 
invariant measure p roughly proportional to 

p{x)ocL^^^{L + x)-^/^{L-xy^/^, x^El, (1.6) 

and that its relaxation time can be bounded from above and below by constant multiples of 
L^'^; see Lemma 2.2 and Lemma 4.6 below for more details. Notice that, although the measure 
p gives uniformly (in L) positive mass to the two attractors x^ = ibL, the drift which pushes 
the random walk away from the saddle x = is proportional to the inverse of the distance from 
the attractors. In particular the naive guess for the mixing time Tmix ~ l/p{0) = 0{L^''^) is 
wrong. 

This heuristics is turned into the rigorous bound (1.4) by using a suitable test function in the 
variational principle that characterizes the spectral gap; see [7, Section 6]. However, it seems 
very hard to give a rigorous upper bound on Tj-d of the same order of magnitude. We obtain a 
bound that can be off by at most two powers of L. 



Theorem 1.1. For any A < 1, 



logTmix / 5 

hm sup — -— ^ - + 2 , 

L^oo log L 2 



The main tool for the proof of Theorem 1.1 is the analysis of an effective dynamics for the 
crossings of the polymer. To describe this, we introduce the variable a G {— 1,+1}'^^, where 
Ol = {—L, . . . , L}\ El denotes the sites with the same parity of L + 1. If r/ is a configuration of 
the polymer, then ry^; / at any x G Ol, and we define a{r]) by ax = sign{r]x)- The projection 
of TT on 5 = {-1, +1}'^^ is then 

,.{a)= Y. vr(r,), (1.7) 

where the sum is over all configurations rj compatible with the signs a. The field i' has non trivial 
long range correlations. Consider the heat bath dynamics for the variables a: sites x G Ol are 
equipped with independent rate 1 Poisson clocks; when site x rings we replace ax by a'^ where 
the new sign a'^ is distributed according to the conditional probability h'{-\ay, y / x), i.e. 
the probability (1.7) conditioned on the value of cr^y , y ^ x. Denote by T^'^j the corresponding 
relaxation time. For this process, the exponent 5/2 can be shown to be optimal. 

Theorem 1.2. For any A < 1, 

L->oo log L 2 
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The proof of Theorem 1.1 and Theorem 1.2 combines several different tools which play a 
prominent role in the analysis of convergence to equilibrium of Markov chains: decomposition 
methods, spectral gap analysis, comparison inequalities, and coupling estimates. An outline of 
the main steps of the proof is given at the beginning of Section 4. 

1.3. Metastability. Recall the definition (1.5) of the two sets Q^, and define the associated 
phases as the restricted equilibrium measures vr^ := 7r(- | $1^), so that (cf. Section 2) 

7r-^(7r++7r-)||=o(l). (1.8) 

In the thermodynamic limit, we expect relaxation to equilibrium within each phase to occur on 
time scales T^^^ such that Tj-ei ^ T^.^^, while on a time scale proportional to T^^i one should see 
the system jump from one phase to the other according to i.i.d. exponentially distributed times. 
A strong indication of this metastable behavior comes from the following theorem. Below, we 
use r/(t) to denote the state at time t of the Markov chain with generator (1.1). The notation 
o(l) refers to asymptotics as L — t- cxd. 

Theorem 1.3. There exists a set S^ C il^ such that tt{S^) = 1/2 + o(l), and that uniformly 
in rj £ S~^ and uniformly in t > 0: 

r'l (r- >t)= e-*/(2^-i) + o(l) , 

where 

T- = inf{t ^ : ri{t) G S'} , S~ = -5+ . 

Remark 1.4. From the proof of Theorem 1.3 it will be clear that the set S*^ is increasing 
w.r.t. the natural partial order among polymer configurations defined in Section 2.5, so that in 
particular the maximal configuration (in the sequel denoted by A) is in S^. 

Here P'' stands for the law of the process with initial state r] £ Vt. By symmetry, Theorem 
1.3 also implies that uniformly in 77 G S~ = —S"^ and uniformly in t > 0, 

t/(2T,,i) 



P^ (r+ >t)= e-*/(2T,,,) ^ ^^^^ ^ 

where t+ = inf{t ^ : r/(t) G 5+}. 
If we define the renormalized process 

ujs = 1 HsTrcx) G f^+) - 1 (r/(sTrci) G VL~) , 

we expect that, starting from any configuration in 17 , {ujg, s ^ 0} converges to the simple 
two-state Markov chain with switching rate (from ±1 to =pl) equal to 1/2, whose spectral gap 
equals 1. Such a strong uniform result seems very hard to obtain for our model. The difficulty 
is that, in contrast with familiar metastability results [18], here there is no dramatic separation 
of time scales: while (1.4) and Proposition 2.6 below imply Tj-d ^ T^^, the ratio T^e\/T^i is only 
polynomially large in L. However, we do have a detailed description of the renormalized process 
when the initial condition is the maximal configuration. Namely, define the maximal element of 
Jl as ?7max = A, i.e. A^^ = x + L for x ^ and A3; = L — a; for x ^ 0, and let T^-^^{e) denote the 
first time t such that ||-Pt(A, •) — 7r|| ^ e. 



Theorem 1.5. For any 6 > 0, uniformly in t ^ L^+ 



^ T^^ H ^ 71"" 



0(1) (1.9) 
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and uniformly in t >0 

_±_ r 1 -I- p — V^rcl 1 _ p — t/Trel -] 

= 0(1). (1.10) 



- 1 _j_ Q — t/Ticl 1 _ g — t/Trel 

- TT^ H TT' 



2 

Moreover, for any e G (0, 1/2) one has 

1 



^^Je)=T,ellog(^-J(l + o(l)). (1.11) 

Remark 1.6. Theorem 1.5 shows in particular that, when the dynamics is started from either 
iT~^ or A, there is no ctii-oj(f phenomenon [13], i.e. the variation distance from equihbrium does 
not fall abruptly to zero, but rather does so smoothly (on the timescale Trei)- That is another 
signature of the metastable behavior of our system and it is in contrast with what one expects 
for the neutral or attractive case A > 1. 

One of the key features of metastability is that, once the system decides to jump from e.g. S""*" 
to S~ , then it does so very quickly on the time scale of the mixing time. We verify that indeed 
this is the case for most starting configurations inside S~^ U S~ . 

Let T denote the random time spent outside S~^ U S~ up to the hitting time of S~: 



7":= / l|^(5)g(5+us'-)c}ds. (1.12) 

Jo 

From Theorems 1.3 and 1.5 one easily deduces that, for most initial conditions in 5^, r^ S> T: 

Corollary 1.7. There exists a subset S^ of the set 5"^ of Theorem 1.3 satisfying tt{S^) = 
l/2 + o(l) such that, uniformly on rj £ S^ , T = o{t^) in probability, i.e. there exists a sequence 
6l tending to zero as L — )■ oo such that for every rj G S'^ , 

r"" [T > Slt-] < 5l. (1.13) 

Along the same lines of the proof of the Corollary, one can establish the weak convergence of 
the renormalized process ujs to the two-state Markov chain, provided that the initial configuration 
is inside a suitable subset of S^ U S~ with almost full measure. We decided to omit details for 
shortness. 

1.4. Organization of the paper. The rest of the paper consists of three sections. Section 2 
starts with standard material and then proceeds with the introduction of some essential tools 
to be used in the proof of the main results, including general results for monotone systems 
that can be of independent interest. This section contains also some new results concerning the 
relaxation within one phase and the properties of the principal eigenfunction of the generator. 
The metastability results are discussed in Section 3. Here, we start with the proof of Theorem 1.5. 
In later subsections we develop the construction needed for the proof of Theorem 1.3. Finally, 
Section 4 proves Theorem 1.1 and Theorem 1.2. This section is broken into several subsections 
corresponding to the various steps of the proof. A high level description of the arguments 
involved is given at the beginning of the section. Finally, 

Notational conventions. Whenever we write o{LP) or 0{LP) for some p € M it is understood 
that this refers to the thermodynamic limit L — )■ oo. Also, we use the notation f{L) = ^{U') 
when there exists a constant c > such that f{L) '^ cU' for all sufficiently large L. For positive 
functions /, ^r, we use the notation f{V) » g{T) whenever liminf2,_!.oo /(-^)/5'(-^) = +oo, and 
/(L) ^ g{L) when limi_s.oo fiL)/g{L) = 1. Also, we write / x (7 if there exists some constant 
c > such that c~^g ^ / ^ eg. 
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2. Some tools 



We begin with some generalities about reversible Markov chains. Then, we recall the definition 
of the polymer dynamics and derive some consequences of monotonicity. Next, we give some 
estimates on convergence to equilibrium in the "plus" phase. Finally, we characterize in detail 
an eigenfunction of L with eigenvalue —gap. 

2.1. Preliminaries. We will consider reversible continuous time Markov chains with finite state 
space X, defined by the infinitesimal generator L acting on functions / : X i— t- M, 

[£/](x) = j;c(x,y)[/(y)-/(x)], (2.1) 

where c(-, •) is a bounded non negative function on X x X satisfying 7r(x)c(x,y) = 7r(y)c(y,x), 
for a probability measure vr on X. In the applications below, the rates c(x, y) will always be 
such that the Markov chain is irreducible and the reversible invariant measure vr is positive on 
X. We refer e.g. to [1, 13] for more details on reversible Markov chains. 

Let vf = P(t;f G •) denote the law of the state wf of the Markov chain at time t with initial 
condition x G X. We shall investigate the rate of convergence of ff to vr. If the initial condition x 
is distributed according to a probability ^ on X, we write v^ = Ylxi^x A*(^)^f ^^ the distribution 
at time t. As usual, one can associate a semi-group {Pt, t ^ 0} to the generator C in such a 
way that [Ptf]{x) = [e*^/](x) = Y^yex ^t{y)f{y)- We also use the notation Pt{x,y) = ff(y), 
and v^ = fiPf. 

The mixing time of the Markov chain is defined by 



T ■ 

-^ mix ^ , , 

xex 
where 



(e) = inf (t > , max \\uf - 7r|| ^ e] , (2.2) 

I = -7y^\f^{x) - V{X)\ 



ll/i — v\ 

is the total variation distance. We shall use the convention that T^^ix stands for Tjnix(2i)- It is 
well known that with this notation one has 

||z/f -7r|| ^e-LVT'n.ixJ^ (2.3) 

for all t ^ 0, where [aj denotes the integer part of a ^ 0. The spectral gap and the relaxation 
time of the process are defined by 

gap = mm — -— , Tj-d = , (2.4) 

f:X^R Var^(/) gap 

where for / : X i— t- M, 

£{fj) = Y, vr(x)/(x)[-£/](x) = 2 E <^y{x,y)[f{y) " f{^)f (2-5) 

x<^X x,yGX 

is the quadratic form of the generator, a.k.a. the Dirichlet form, while Var,r(/) stands for the 
variance 7r(/^) — 7r(/)^. Thus, gap is the lowest non zero eigenvalue of —C The following bound 
relating total variation distance and relaxation time is an immediate consequence of reversibility 
and Schwarz' inequality: 

||^.f-7rK^e-*/^-'v/Va^, (2.6) 
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where /(o") = ^{a)/'K{a) and /U is a probability on X. Another standard relation between total 
variation and relaxation time is the identity 

gap = — lim - logmax ||ff — vf\\- (2-7) 

t^foo t x,y 

Combining (2.3), (2.7) and (2.6), one can obtain the following well known relations: 

Trc\ ^ Tlnix ^ (l — log7r*)Trci , where tt* = min7r(x) . (2.8) 

' xex 

2.2. A general decomposition bound on the spectral gap. We shall need a continuous 
time version of a general decomposition bound obtained by Jerrum et al. [11]. Consider the 
continuous time reversible Markov chain defined by (2.1). Suppose the space X is partitioned 
in the disjoint union of subspaces Xi, . . . , Xm, for some m £N and define the generators 

[^if]{x) = '^Ci{x,y)[f{y)- f{x)], Ci{x,y) = c{x,y)l{y e Xi) , x e Xi . 
yex 

Then d is the generator of the Markov chain restricted to Xi, its reversible invariant measure 
being given by ttj = 7r(- 1 Xj). Let Amin denote the minimum of the spectral gaps of the Markov 
chains generated hy Ci, i = 1, . . . , m. Next, let C denote the infinitesimal generator defined by 

m 

[Cip]ii) = ^c{i,j)[ip{j) - ip{i)] , 
i=i 
for (p G M™, where 

c{i,j)= ^ TT{x\Xi)c{x,y). 
xGXi , y&Xj 

This defines a continuous time Markov chain on {1, . . . ,m} with reversible invariant measure 
7f(i) = Tr{Xi). Let A denote the gap of this chain. A straightforward adaptation of [11, Theorem 
1] yields the following estimate. 

Proposition 2.1. Define 7 = maxj max^gXj J2y£X\x- ^(^' v) ■ Then, with the notation of (2.4), 

. r A A Amin 1 /^ „^ 

2.3. Killed process and quasi-stationary distribution. Here we recall some standard facts 
about killed processes, their generators and quasi-stationary distributions for reversible Markov 
chains; we refer to [1] for an introduction. Given a reversible Markov chain with generator C as 
above, and a subset T C X, we consider the process killed upon entering T, with sub-probability 
law defined by 

u^'^{B) = F^'ivt eB; Tr>t), x G ^ , (2.10) 

where B C X, vt denotes the state of the Markov chain with generator C at time t, P^ denotes 
the law of the process started at x, and tt denotes the hitting time of the set F. The associated 
semi-group Pf is given by 

[P[f]{x) = [e'^'f]{x)=Y,i.f{y)f{y), xG^, (2.11) 

where the killed generator C^ satisfies, for every x G V^: 

[C^f]{x) = [C{f lr^)]{x) = [Cf]{x) -Y^c{x,y)f{y) . (2.12) 

y&V 



DYNAMICS OF A POLYMER WITH REPULSION 



We assume that Pf is irreducible. Then C is a negative definite, self-adjoint operator in L2('7r), 
and its top eigenvalue — 7r is characterized by 

7r = mm = mm , (2.13) 

f:X^R, 7r(/2) f:X^R, Tt{P) 

/lr=0 /lr=0 

where we use (•, •)^ for the scalar product in L2(7r), and £{f, f) is defined by (2.5). 

Let gr denote the (unique, positive on V^) eigenfunction of £^ associated to — 7r. Extending 
gr to all X G X by setting gr{x) = for x G T, one defines the quasi- stationary distribution u^ , 
i.e. the probability on X given by 

.,(!,) = =M?eM, ,,x. (2.14) 

An equivalent characterization of it is as the limit 

uriB) = lim P^(t;t G S | rr > t) , (2.15) 

where B C X, and the chosen initial point x G T'^ is arbitrary. The fundamental property of 
the quasi-stationary distribution is that, starting from i^r, the hitting time rr is exponentially 
distributed with parameter 7r: 

F''^{TT>t) = e-^^\ (2.16) 

where P^"" stands for the law of the of the process when the initial state is distributed according 
to fr. Another way of expressing quasi-stationarity is y Pf = e"''''" V , for all t ^ 0. 
A general property of 7r (cf. Lemma 3.1 below) is that 7r > gap7r(r). 

2.4. Polymer model. Let O = Q2L stand for the space of all lattice paths defined in the 
introduction. A partial order in il is given by 

rj <^ 7]' ^^ 7],^ i^ 7]'^ , x = -L,...,L. (2.17) 

Given C,^ ^ Q such that C ^ C we define the restricted space il'''^ of all paths rj ^ Q, such that 
C ^ V ^ C- The dynamics is defined by the continuous time Markov chain with state space O''''', 
with infinitesimal generator C^'^ given by (1.1) where the rates rx,±{r]) are replaced by 

^^;i(^) = r.,±(r?)l(??"'^ G n'^^^) (2.18) 

This process is the heat bath dynamics associated to the probability measure tt^^^ on 0^'^ 
defined as in (1.2) with the normalization now given by 

zk'''= E ^"'^''^- (2-19) 

Equivalently, vr2^ = 7^21^' I ^^'^)- This is referred to as the polymer model with top/bottom 
constraints (C is the bottom, ^ is the top). For simplicity, when no confusion arises, we often 
omit the superscripts A, C, ^ and the subscript L from our notation in what follows. We write 
Vf for the state of the Markov chain at time t when the initial configuration is some ij, and let 
u^ denote its distribution. When the initial condition rj is distributed according to a probability 
measure /i on fi we write u^ as in Section 2.1. 

Note that the generator C^'^ can be written in the form (2.1) by setting c{r], rj') = r^'^{r])l{ri' = 

r/^' ), and n = 'n'21 is reversible. While this holds for every value A > of the parameter 
describing the strength of the interaction, we will only consider the case A < 1 below, which 
corresponds to a strictly delocalized regime for the polymer. 
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The minimal path V and maximal path A for the order (2.17) are defined by V^; = —x — L 
for X ^ 0, Va: = -L + X for X ^ 0, and A = -V. Clearly, ii ( = V and ^ = A, then Q'^'^ = Q,. 
This case is referred to as the polymer model with no top/bottom constraint. 

The following well known estimates will be often used in our proofs. We refer e.g. to [9, Section 
2] for the proof of Lemma 2.2 below, as well as for other known properties of the delocalized 
equilibrium measure. Let Z2L = ■^2L denote the partition function (2.19) with no top/bottom 
boundaries and write ZJ^ = Z2I for the partition function (2.19) with ^^ = A and ( given by 
the minimal non negative element of fi, i.e. Cx = if x G El (x has the same parity as L) 
and Cx = 1 if 2; G Ol {x has opposite parity w.r.t. L), i.e. Z^ is the partition function of the 
polymer with a horizontal wall at height zero. Recall that N = N{r]) stands for the number of 
zeros in the path r/ lying strictly between —L and L. Considering reflections of the path between 
consecutive zeros one obtains 

2^2'2' = ^2f- (2-20) 

Lemma 2.2. Consider the polymer with no top/bottom constraint with A € (0, 1). There exist 
constants ci = Ci{\) > 0, z = 1, 2 such that 

2-'^^Z^L ~ ci L^3/2 , (2.21) 

and 

7r(iV(r/) >k)i^C2 e'^l^^ , V A; G N . (2.22) 

An immediate implication of (2.20) and (2.21) is that 

2-2L^+,A ^ ^^ ^-3/2 ^ ^2.23) 

for some constant c+ > as soon as A < 2. Moreover, (2.21) and (2.23) imply the bounds 
vr (r/j/ ^ Vy G {-L, . . . , x} , and ??x = 0) x vr (77^ = 0) 

xL3/2(L + x)-3/2(L-x)-3/2, (2.24) 

for every x G Ei^. 

2.5. Monotonicity. An important property satisfied by the Markov chains introduced above 
is the monotonicity with respect to the partial order (2.17). A convenient way of stating the 
monotonicity property is that there exists a coupling P of the trajectories of the Markov chains 
corresponding to distinct initial conditions such that if r/ ^ r/' then P almost surely v^ ^ vl for 
all t ^ 0. More generally, one can define a coupling P of trajectories corresponding to distinct 
top/bottom constraints and distinct initial conditions such that if C ^ C'l C ^ ?') ^md 77 ^ ??', 
then P almost surely v^ ^ Vf for all t ^ 0. Recall that a function / : $7 1— t- M is said to be 

increasing if /(??) ^ fiv') whenever rj ^ rj' . An event A is increasing if the indicator function 1a 
is increasing. The monotonicity property of the dynamics implies the so-called FKG property 
of the equilibrium measures n = t^2L ■ ^^'^ every pair of increasing functions /, 51 : ri 1— )■ M, one 
has 7r{fg) ^ Tr{f)Tr{g). We refer to [7, Section 2] for a more detailed discussion of the monotone 
coupling and the consequences of monotonicity. 

Lemma 2.3. Let fi be a probability on Q, and write f{rj) = fj,(ri)/TT(ri), and ft{ri) = f^(r/)/7r(r/), 
t > 0. If f is increasing then, for every t > 0, ft is increasing. As a consequence, there exists 
an increasing event A such that 

||^r-7r||=^f(^)-7r(A). (2.25) 
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Proof. Write ff(r/) = X^„„gQ M(%)-Pt(%5 ^)) where Pt{-,-) stands for the kernel of the Markov 
chain. Reversibility then gives 

Mr]) = Y^ /(r?o)vr(r/o)Pt(%,r/)/7r(r/) = J] f{r]o)Ptiv,Vo) ■ (2.26) 

Next, let P denote the monotone coupling introduced above and let E denote expectation w.r.t. 
P. Then, (2.26) coincides with E[/(fJ')], and if r? ^ f]' , 

fM) - fM = nfivf) - /(^I')] = E[/(^;'') - /(^I') ; vf ^ t;?] . 

Thus, ft is increasing whenever / is. Finally, it is well known that the total variation distance 
can be written in the form (2.25) where A = {r] : i^f{A) ^ tt{A)}. Since A = {ft ^ 1}, A is 
increasing whenever /is. D 

Lemma 2.4 compares arbitrary initial conditions to the extremal initial conditions. Lemma 2.5 
states a useful sub-multiplicativity property satisfied by extremal evolutions. For lightness of 
notation, we state these results only in the case of no top/bottom boundaries, i.e. C = V, ^ = A, 
but the same applies for general C, ^ with exactly the same proof. 

Lemma 2.4. For any t > and any rj, r( G ^: 

II n »?' 1 1 ^ /I r 2 1 1 A V 1 1 

If* ~ ^t II ^ 4/v iiz/j — vi II . 

As a consequence, 

gap = - lim - log \\v'^ - u'^ II . 

i— >oo t 

Proof. Let P denote the monotone coupling as above. Then, 

L-l L-1 

x=-L+lh=-L 
^4L2||iyf -Z//II. 

The second point follows from the first one and the classical characterization (2.7) of the spectral 
gap. D 

Lemma 2.5. For any s,t^ 0, 

||,yA _ , .V II < ||,,A _ Vll II A _ Vll 

Irt+s '^t+sW ^ W'^t '^t II IFs '^s II • 
Proof. With the same argument of Lemma 2.3, for some increasing event A 

H^s-i^]!+s\\ = <s{A)-v^t+s{A). 
Let p be a coupling beween v^- and u^ at fixed time t ^ 0. Then 

v^^M) - ^ts{A) = [{u^iA) - u'l{A))dp{n,a) 



(i.,^(A)-<(A))l(a/77)dp(r/,a) 

^(^^(A)-^r(^M^/^) 

To conclude, we take p as the maximal coupling, i.e. such that p{cT ^ rj) = \\v^ — u^ \\. D 
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2.6. Relaxation in one phase. Here we obtain results concerning tlie polymer dynamics in 
the phase tt^ defined after (1.5); cf. Proposition 2.6 below. Then, we show that the polymer 
started at the maximal configuration A relaxes first to the restricted equilibrium vr"'' in a time 
0(L^^ ) for arbitrarily small 6 > 0, while it takes much longer to reach the full equilibrium vr; 
cf. Lemma 2.7 and Lemma 2.8 below. 

Recall the definition (1.5) of the subspace i7^ C ^, where L ^ i, and i diverges as L — ;■ oo; 
see (2.29) below. The corresponding restricted equilibrium is given by vr^ = Tr{-\Q~^). Note that 
this is a particular instance of the polymer equilibrium vTg); with top/bottom boundaries: the 
top is ^ = A while the bottom (" = C(^^) is the lowest element of ri"*". Similarly, one defines 
il^ = — r2+, i.e. use (1.5) with 7]^ > replaced by r/^; < 0, and the equilibrium vr" is defined 
accordingly. 

Since A < 1, the equilibrium bounds (2.24) imply that 

7Tin+) = 7r{n~) = ^ + oir^/^); (2.27) 

see e.g. [7, Section 2]. In particular, if i diverges as L — ;■ oo, then 

llvr- ^(7r++7r-)|| =o(l). (2.28) 

What follows depends only marginally on the precise dependence of ^ on L, provided that 
L ^ i ^ 1. For the sake of simplicity we shall fix its value as 

^(L) = (logL)i. (2.29) 

This choice turns out to be convenient in the proof of Proposition 2.6 below, but we point out 
that any choice of the form i{L) = 0{L^) for small e > would be sufficient to obtain the same 
conclusion with a little more work. 

We start by establishing a mixing time upper bound for the dynamics constrained to stay in 
0+, i.e. the process evolving with bottom boundary given by C = C(^^)- To avoid confusion we 
shall write /ij' (instead of f^) for the law at time t of this Markov chain with state space il^ and 
initial condition rj. We write C~^ for its generator and gap+ for the associated spectral gap. 

Proposition 2.6. For every e > 0, there exists Lq = Lo(e) such that for all L ^ Lq, for all 

t ^ and all initial conditions rj £ Q~^ : 

Wfil - TT+IK 4L2 exp (- t/L^+') . (2.30) 

In particular, gap+ ^ L^^^^ . 

Proof. The last statement follows from (2.30) and (2.7). To prove (2.30) we establish that for 
every e > 0, there is a constant Lq = Lo(e) > such that, taking T = L^+'^, we have 

||;U^-^^|| ^ 1-L"S (2.31) 

for all L ^ Lo(e), where (" stands for the minimal element C = C(^^) of ^^- Once (2.31) is 
available, we obtain (2.30) (with a new value of e) from Lemma 2.4, since Lemma 2.5 (which is 
also valid for the restricted dynamic) and (2.31) imply 

\\u^ //IK (\\u^ u<\\V'^^^ 

^ exp {-[t/T\L-')) ^ 2 exp {-t/L^+^') , 

for any L large enough. 

To prove (2.31), we divide the sites x G {—L-, ■ • ■ -.L} in three overlapping regions: 

h = {-L, . . . ,-L + f) , h = {-L + I, . . . ,L - I] , s.nd h = {L - e , . . . ,L] , 
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where I = i{L) is given by (2.29). Let T2 = L^+'^i, Ti = L^^ with some ei > such that 
T = L2+^ ^ T' := 2T2 +Ti. We shall prove (2.31) with T replaced by T' (this imphes the 
claim since by Lemma 2.4 the left hand side of (2.31) is monotone as a function of T). Call 
/i^°''^ the law of the "censored" process obtained as follows. Start from r/o at time and, for 
time t £ [0,72] reject all the updates involving x ^ l2- For time t £ (T2,T2 + Ti] reject all 
updates involving x ^ IiU I3, and for time t G (T2 + Ti, 2T2 + Ti = T'] reject all the updates 
involving x ^ l2- From the Peres- Winkler censoring inequality [19, Theorem 16.5] one has that 
/ij,f stochastically dominates fXj,,. Similarly, fijif is stochastically dominated by ^u^,. On the 
other hand, as in Lemma 2.3, one has 

where A C Jl"'" is an increasing event. Therefore, 

and the lemma follows once we show that 

ll/i^f - ;u^f II ^ 1 - L-^ (2.32) 

To prove (2.32) we shall couple the two configurations r]rpf, rjj^^ with law /U^f , fiqlf respectively. 
From the analysis of the polymer with the wall [7, Section 4] , it is not hard to infer that uniformly 
in the boundary values at —L + i — 1 and L — i + 1 the system evolving in the region I2 has 
a mixing time 0(L^ log L). Therefore, after a time T2 = L^+^i, up to 0{L^^) corrections for a 
large constant p > (for L ^ Lo{p)), for any event E, firp' (E) coincides with the equilibrium 
probability of E in I2 with boundary conditions A at —L + i— 1 and L — i + 1. The same 
applies to iiY^{E) provided the equilibrium is taken with boundary conditions C, at —L + i — 1 
and L — i + 1. Choose the event E that the configuration is minimal (in Q'^) at both —L + i'^ — 1 
and L — i'^ + 1 (i.e. r/_2,+£2_i = f/L_£2^i = if L — ^^ is odd, and ??-__l+£2_i = f/L_£2^i = 1 if 
L — l"^ is even). From known equilibrium estimates [9, Section 2], it is not difficult to show that 
at equilibrium, with either of the two boundary conditions considered above, the probability of 
E is bounded below by ci i~^ for some constant ci > depending only on A. Therefore, using 
an independent coupling in the time-lag [0, T2], we have that the event E occurs for both rjrp' , 

r/^''^ with probability at least c£~^'^. 

Next, conditioned on the event E we see that from time T2 up to time T2 + Ti the two 
processes evolve (in the regions Ii and I3 only) with the same boundary conditions (equal to the 
minimal configuration at both —L + i'^ and L — i"^). Since Ti = L^ and the mixing time of the 
system in the regions /i and Is (which evolve independently) is certainly at most e '^ •*, with 
probability 1 + 0{L~p) (conditionally on event E) we have that the two configurations coincide 
in both regions /i , I2 at time T2 + Ti. Therefore if we let the system run only in I2 now for an 
additional time T2 we have a probability close to 1 to have the two configurations coinciding 
everywhere. It follows that, conditionally on the event E, there is a probability of, say, at least 
1/2 of no discrepancy between r/^f and r/^f . Therefore, letting P denote the coupling described 
above, 

II A.c C,c\\ ^ TTT)/ A.c _i_ C,c\ 

ll/n^', - n)^, II ^ F{r]j,] yt: rjr^^ ) 

^ P(??^;" / 4f I E)FiE) + 1 - FiE) 

^i-\f{e). 

Since F{E) ^ cjt'^ ^ 2/L^ for L ^ Lo(e), this implies (2.32). D 
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We now go back to the model with no top/bottom boundaries, that is the law v^ corresponds 
to the evolution with ^ = A, C = V. The next result is crucially based on estimates obtained in 
[7, Section 6] for the delocalized regime A < 1. 

Lemma 2.7. Uniformly in t ^ L^/^(logL)^^, 

u,^{n+) = 1 + 0(1). 

Proof. Define the event A = { "^^=-1^^ ^ L^'^(logL)^^ }. Proposition 6.2 in [7] proves that 
t'^(^) = o(l) uniformly in t ^ L^/^(logL)~^. Since Q~ C ^1 we have 

u^{n-) = o{l), t^L^/\logL)-\ (2.33) 

Next, let us check that 

iy^{n+\{n-y)^TT{n-^\{n-y), t^o. (2.34) 

To this end, observe that since {^^^ is increasing, using Lemma 2.3 the function 

is increasing. Since Q^ C i^^)^ is increasing, with the FKG property for n, this implies (2.34). 
From (2.33) and (2.34) we obtain 

iy^{n+) ^ (1 + o(i))z>t^(^+ 1 in~y) ^ (i + o(i))7r(f)+ 1 in-y) = i + o(i) , 

where the last bound follows from (2.27). D 

Lemma 2.8. For any e > 0, uniformly in t £ [L^^''^, L^'^(logL)~^] : 

||z.,^-7r+||=o(l), 
where vr^ is defined by vr^ = it{- \ Q'^). 
Proof. Using Lemma 2.7 it is enough to prove 

Hi.\n+)-n+\\ = oii), 

uniformly in t G [L^'^^,L^/^(logL)~^]. Consider the function / : fi i— )■ M given by 

•^("^ = '^"^") vr+(a) = '^"^"^ ^(oT 14W) ■ 

Since fi^ is increasing. Lemma 2.3 shows that / is increasing. Therefore, the event 

A = {aen+ : z/f (a | n+) > Tr+{a)} 
is increasing. Using monotonicity we have 

u^{A 1 0+) = ^f (^)/^f (f^+) ^ f,nA)/u^{n+) , 

where n^ denotes the evolution constrained to stay in il"^; see Proposition 2.6. Therefore, 

Hi-\n^)-^^ = -nA\n^)-n^A)^J^^-n^{A). 

The conclusion now follows from Lemma 2.7 and Proposition 2.6. D 

The full power of Lemma 2.8 will be seen in the next sections. One of its consequences is the 
fact that the mixing time Tmix can be bounded in terms of the relaxation time via 

Tmix^i^+" + cT,eilogL, (2.35) 
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for some constant c > 0. Indeed, (2.35) follows quite easily from Lemma 2.4, Lemma 2.8 and 
(2.6); see Lemma 4.1 below for a more subtle application of the same reasoning. Note that, since 
Trci S> L^/^""^, the bound (2.35) improves considerably the standard estimate (2.8) by replacing 
the factor — logvr* = 0{L) with a factor O(logL). 

2.7. Characterization of the principal eigenfunction. What follows refers to the model 
with no top/bottom boundaries. Recall that A, (resp. V) denotes the maximal (resp. minimal) 
configuration in Q. A function (/ : fi i— )■ M is called antisymmetric ii g{—rj) = —g{r]) for all rj £ il.. 
The following result gives a precise characterization of one eigenfunction corresponding to —gap. 

Proposition 2.9. There exists an increasing antisymmetric eigenfunction g of C, such that 
WaWhiiTT) = 1- It satisfies 

C-a = -gap g- (2.36) 

Moreover, when L tends to infinity 

5(A) = |bl|L^=l+o(l). 
and 

||5-(W-ln-)lkw = o(l). (2.37) 

Proof. By decomposing l/\ — ly on a basis of eigenfunctions of C one sees that 

,. Pt{l^ - Iv) 
g := lim - — — . 

t->oo ||Pt(lA - lv)||L2(7r) 

is an eigenfunction with unit L2 norm. It is increasing and antisymmetric as -Pt (1a — Iv) is 
antisymmetric and increasing for all t (Pf preserves monotonicity and symmetries). To prove 
(2.36), it suffices to show that the projection of 1a — Iv on the eigenspace of C associated to 
—gap is non-zero. To do so, first observe that by reversibility, 

^||Pt(lA-lv)||LiW=2||l^r-'^t''l|- 

Then, by the second point of Lemma 2.4 

lim -log||Pt(lA - lv)||L2(7r) = lim - log ||Pt(lA - lv)||Li(7r) 

= hm - log II ff - v^ II = -gap , 

r^oo t 

(where we used equivalence of the norms in finite dimensional spaces) . 

We now estimate the Lqo norm of g. Let e > be small and fixed, and let to be such that 
L^"*"^ < to < L^''^~^- The function g is an eigenfunction for P^,, = exp(to^), with eigenvalue 
g-togap^ Therefore 



vr 



e-*os-P5(A) = E [5«)] ^ 7r+(<7) + 25(A)||< 

where the last inequality follows from the fact that for any two measures ^, v and any function 

/, 

|M/)-K/)I^2||/||l^||^-i.||. 

Hence, by Lemma 2.8 and the fact that gap~^ ^ to (cf- (1-4)): 

5(A) ^ _,,g,p ""^jf I ^ = vr+(5)(l + 0(1)). (2.38) 

1 1 TQ II 
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Moreover, by symmetry and Jensen's inequality 

so that 

Therefore g{A) ^ 1 + o(l) (and it is trivial to notice that Hs'IIloo ^ ll5'llL2{7r) = !)• 

We turn to the proof of (2.37). First notice that by (2.38) one has 7r+ (5) ^ (l+o(l))fif(A) ^ 1+ 
0(1) so that 

^+(5) = 1 + 0(1). (2.39) 

Next, we prove that the variation of g within Q, is small. Let C^ be the generator of the 
Markov chain restricted to fi"*", as in Proposition 2.6. The associated Dirichlet form £~^ is, for 
a function / G L2(vr+), 

£^if,f) = l E ^^(V)civ,v')[f{r]')-f{r])?, 

where c{ri,r]'), rj,ri' G i7+, denote the transition rates, which coincide with those of C If L is 
large enough, gap ^ L^^/'^+^ by (1.4), while gap"*" ^ L~^~^ by Proposition 2.6, so that 



L-5/2+- ^ gap = £{g) ^ 7r{n+)£+{g\n+) 

^ gap+7r(il+)Var^+(5(|Q+) ^ -L"^"^ Var^+(gf|!^+) , 

where we let g\ ^ := glQ+. Therefore, one has 

Var,+ (<7|n+)^4L-i/2+2^ (2.40) 

and the same is true for Yav^- (giQ-) , by antisymmetry. 
Next, 

115 - ln+ + lf^-||Li(7r) ^ l|ln+(5 - l)llLi(7r) + W'^Q- id + l)llLi(7r) + IblnVC^+UH-) IIliCtt) • 

The first two terms of the right-hand side are equal by symmetry. Adding and subtracting 
7r~^{g), and using Schwarz' inequality. 



l^n+ig - 1)IIliW ^ <^ ) 1^ (5) - 1| + VVar,+ (5| 



ofl) 



where the conclusion follows from (2.39) and (2.40). The third term ||5ln\(r2+uc-)llLi{7r) is 
smaller than ||5'||loo^(^ \ (^^ U n^)) = o(l). □ 



3. Metastability 

In this section we first prove Theorem 1.5, which is mainly a consequence of the technical 
lemmas of the previous section and then move to the proof of Theorem 1.3 and its corollary. 
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3.1. Proof of Theorem 1.5. We use the notation T = Li^^^ . Equation (1.11) is an easy 
consequence of (1.9). Indeed, assuming (1.9), for t >T, one has 



\Pt(A, 



TT\ 



r 1 + e"*/^--<=i 



7r+ + 



1 



3 — t/Trcl 



-vr 



1 



7r|| + o(l) = -e-*/^-i + o(l). (3.1) 



To prove the rest of the result, one first shows that proving (1.9) reduces to prove (1.10): 

p '•/ -^ rcl -| 

2 J 



r 1 + e^*/^--<=i , 1 

^ T^ H 



< WrPt 



t-T 



vr+Pt-T + 



vr+Pt- 



T 



1 + e~*/^- 



■vr+ + 



1 



(3.2) 



The inequahty is just triangular inequality, combined with the observation that u/pPT-t = '^t'- 
The first term on the right hand side is smaller than ||z/^ — 7r^|| (as Pt^i contracts the norm) 
which is itself small, by Lemma 2.8 and the definition of T. It remains to estimate the second 
term i.e. to prove (1.10). 

To do this, we use the fact that the density of vr"*' w.r.t. vr is very close to g, the eigenfunction 
described in Proposition 2.9, so that the density of 7r"'"Pt must be close to Ptg. Using reversibility, 
one can express the densities as follows: '^^ * = Pt^:^- Then we rewrite the second term in 
(3.2) as an Li norm (omitting a harmless factor 1/2) 



Pt-T- 



d7r+ 



1 



dvr 27r(f]+) 



Lf7+ 



+ lf 



27^(J^^ 



-t/T-cal 



< 



Pt- 



d7r^ 



1 



T 



dvr 27r(^+) 



Ln+ 



+ lc 



1^+ - lf7 
1 



Li(7r) 



27^(J7^ 



g-t/Trcl„ 



1 



27^(fi^ 



3 — t/Trol 



Ln+ 



Li{7r) 



The last term above is small by Proposition 2.9. From (2.27) we know that 2tt{VL'^ 
One can then estimate the first term 



9) IIliCtt)- 
= 1 + 0(1). 



Pt- 



d7r^ 



1 



dvr 27^((7^ 



(if7+ + in-; 



1 



e-t/T.ci^ 



2t:{VL 

< II (1^+ + 1^-) - l||Li(7r) + 



Li(7r) 



Pt 



d7r+ 
dvr 



Li(7r) 



+ 0(1) 



(3.3) 



where we used the triangular inequality, the fact that Pfl = 1, and 



Pt-TQ 



e-i^-T)IT..,g 



^ — t/Trcl 



+ o{l))g.. 



which follows from T = o{T^^\). On the right hand-side of (3.3), the first term is small by (2.28) 



and the second is bounded by 



| d7r+ 
dvr 



1 — fi'||Li(7r)) which is small by Proposition 2.9. 



D 



3.2. Proof of Theorem 1.3. Theorem 1.5 gives some intuition on why the result should be 
true, and it will be used to determine the time of the jump from one state to the other. However, 
one needs another key ingredient to get the result, namely the description of the quasi-stationary 
distribution. The reason for this is that starting from the quasi-stationary distribution, a killed 
process dies exactly at exponential rate; see Section 2.3. Therefore, most of our efi'ort will focus 
on stochastic comparison with quasi-stationary distribution. Let us first give a brief roadmap 
to help the reader through the proof of Theorem 1.3. 
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Step 1. The sets S of Theorem 1.3 for which we have the desired exponential hitting time 
description are constructed by successively refining a first attempt. One first defines S^' as the 
sets of polymer configurations where the eigenfunction g in Proposition 2.9 is positive (negative) 
and one verifies that their equilibrium probability is ^ + o(l). Then one examines the Dirichlet 
problem associated to the process killed in S*"'" (5*^'+) and one proves that the corresponding 
eigenvalue 70 is of the same order as the spectral gap apart from a crucial unspecified multiplica- 
tive factor in [1/2, 1]. Similarly one verifies that the corresponding quasi-stationary measure is 
very close to vr^, the equilibrium measure vr conditioned to be in S^'^ {S^~). In this way we 
get the exponentiality of the hitting time of e.g. S*"'^ starting from 7r+ with a rate which is, 
modulo a multiplicative factor in [1/2, 1], the spectral gap (see Lemma 3.3). 

Step 2. Next one appropriately defines new sets S^' C S^' in order to guarantee that this time 
the corresponding Dirichlet eigenvalue 71 is equal to (2 +o(l))gap, and that the hitting time of 
S*^'^ starting from equilibrium conditioned to S^' is exponential (with the correct rate). Again 
one of the key points is to show that vr"^ is close to the quasi-stationary distribution associated 
the process killed on entering 5"^'^, and that the equilibrium probability of S*^' is still ^ + o{l). 

Step 3. Finally, one defines the final sets S*^' C S"^' in such a way that: a) the hitting time of 
S*^'^ starting from any configuration in S*^'^ (and not just from the conditional equilibrium) is 
also exponential with the correct rate ^gap; b) the equilibrium probability of 5^'^ is still ^+o(l). 

It is now time to begin the implementation of the above strategy. Let 

50.+ := {r^ G f), 5(^)>0}, (3.4) 

where g is the eigenfunction defined by Proposition 2.9 and S^'^ = —S^'^. From Proposition 2.9, 

III50.+ - ln+llLi(7r) = 0(1)- (3.5) 

In particular, 7r(5°'^) = 1/2 -|- o(l). Let S~ C S^~ be a decreasing event. We consider the 
quasi-stationary distribution v^ := ug- of the process killed when it hits S~ . Let P^* = P^ , 
resp. C* = C^ , denote the semi-group, resp. the generator, associated to this process (see 
Section 2.3), —75- = —7 be the largest eigenvalue of C* and r" = Tg- be the hitting time of 
S-. From (2.16): 

F""^ {t~ > t) = e-^K (3.6) 

Our first step is to prove that if S^ has non-negligible measure, then 7 is of the same order of 
the gap. More precisely: 

Lemma 3.1. For any S~ C S^~, one has it{S~) ^ 7Trci ^ 1. 

Proof. The bound tt{S~) ^ 7Trci is rather standard, but we include its proof for the sake of 
completeness. Let /o = gs- denote the minimizer in the variational principle defining 7 = 75-; 
see (2.13). Then 

Var„(/o) = (/o,/o).-</o,l(5-)=>' ^ (/o, /o). vr(5^) , 

where we used the Cauchy-Schwarz inequality for (|/o, l/g-y'j . Therefore 

Sifo,fo) . fc-^Sifo,fo) ^ /c'-^ 
7r(/o) Var^(/o) 
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As for the bound 7Trei ^ 1, 75- being a non-decreasing function of S~ (for the inclusion), 
it is sufficient to prove the result for the maximal case S~ = S^'~ . Let g be the eigenfunction 
defined in Proposition 2.9. From (2.12), for all r] G {S^'~Y 

-i^*9\^go,-)c)iv) = -{^9){'n) + ^ c{T],r]')g{r]') ^ - iCg)i7]) = gap g{r]) , 

where we use the fact that g{r]') < for t]' G S*"'^. Plugging this into (2.13), and using 
ff|(so,-)c ^ 0' one gets 



D 



7 ^ ]-2 N ^ gap . 



Next, we prove that the quasi-stationary distribution z^"'' for the process killed on S is very 
close to 7r+ if S~ has probability close to 1/2. 

Lemma 3.2. Uniformly for all decreasing events S~ C S^'~ , 

Ik^-^^ll ^ (2-47r(5~)) + o(l). 

Proof. We use triangular inequality to get 

IIZ.+ - vr+K lk+ - <■ \iS-n\ + ||vr+ - n{- |(5-)'=)||. (3.7) 

We start with the first term. First, from (2.15) one has the characterization 

u = lim . 

Since the operator P^ preserves monotonicity {S~ is decreasing), arguing as in [19, Lemma 16.6], 
the density ^^' is seen to be an increasing function for every fixed t ^ 0. Hence, passing to 
the limit t — )■ 00, dz^^/dyr is an increasing function. Therefore, 

A := I, G iS-y, such that ^^^^^f^'^^ > 1 

is an increasing event. From standard properties of the total variation distance 

\\u+-7ri-\iS~r\\ = ,.+ iA)-niA\iS~r). 

We shall prove that i^~^{A) is smaller than 7r~^(A) + o(l) by the use of monotonicity and a chain 
of comparisons. Recall the notation T = L^"*" [5 G (0,1/4)). We first compare v'^ to u^ Pt'- 
remark that 

v+Pj, = u+p* + z/+(Pr - Pt) 
where the two terms of the decomposition are positive measures. From quasi-stationarity one 
has u'^P^ = e~'^'^u^ and therefore the total mass of the second term above is 1 — e~^'^ . Hence 



1 d|,^:^(ft-pf)i _ dj^ 

" dTT dvr 



^ l-e-T^ = o(l). (3.8) 

Li(7r) 



" 2 

The last equality comes from Lemma 3.1 and the fact that Tj-ci ^ T. Next, from Lemma 2.3, 
6;^Pt stochastically dominates v^Pt so that [i/^Pr](yl) ^ [5/\Pr](yl). Hence, from Lemma 2.8 
and (3.8): 

u+{A) ^ u+Pt{A) + 0(1) < 5„Pt{A) + 0(1) 

^ TT+{A) + \\5;,Pt - vr+ll + 0(1) = 7r+(A) + o(l). 
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Therefore, going back to (3.7) 

\\v+ - TT+II ^ v+{A) - ^{A \{S-f) + IK+ - 7r(- \{S-r)\\ (3.9) 

^ vr+(A) - ^{A \{S-f) + IK+ - 7r(. |(5-)^)|| + o(l) ^ 2||vr+ - vr(. \{S'r)\\ + o(l) . 

To estimate the right-hand side of (3.9), notice that 

|K+-7r(. \{S-r)\\ ^ ||7r+-vr(. |(50'-)^)|| + M- \{S''-f)-^{- |(5-)^)||, (3.10) 

and the first term is o(l) by Proposition 2.9. Moreover, since S^ C S^~ 

||,r(. 1(5°'-)^) -7r(. |(5-)^)|| = ^('5")') -^((-5°'-)^) 



l/2-7r(5-) + o(l) 



^_^(^_) ^l-2vr(5-) + o(l). (3.11) 

Combining (3.9), (3.10) and (3.11), the desired result follows. D 

Now one uses the fact that v^ and vr^" are close in total variation distance to estimate the 
jumping time to S~ starting from either A or from vr"'". For the rest of this section, one defines, 
in analogy with r", the hitting times r*'~ (resp. r*'+), {i = 0, 1,2) of the sets 5"*'" (resp. 5*'"'') 
to be defined. 

Lemma 3.3. Uniformly for all t^ Q and all decreasing S^ C S^'^ , setting 7 = 75- : 

(i) |P^ [t~ >t]- e-T*| < (2 - 47r(5-)) + o(l) 
(ii) |P^+ [t- > t] - e-T*| < (2 - 47r(5-)) + o(l) 

In particular, for S~ = S^'~ , setting 70 = 750,- .' 

(iii) P^ [t°'^ > *] = e^^"* + 0(1) 
(iv) P'^^ [r° - > t] = e"T"* + 0(1) 

Proof. Item (ii) follows from (3.6) and Lemma 3.2. Indeed, 

IP'^^ [t- >t]- W"^ [t- > t] K lk+ -v+\\- 

For item (i) (lower bound), we use the fact that S^ is a decreasing event to get that 

P"" [t- >t] ^ F"^ [t- >t]= e"^*. 

For the upper bound it is sufficient to prove the result for t ^ T = L^+ , since 7^^ S> T (a 
consequence of Lemma 3.1 and (1-4)). One defines 

r~ = inf{t ^ T, 7]{t) eS~}. 

Then, by the Markov property and quasi-stationarity 

P^ [r- >t] < P^ [r~ >t]= F^'^^T [t- >t-T] ^ \\S;,Pt -i^^\\+ e'^^*"^) 

^ \\6aPt - vr+ll + ||7r+ -iy+\\+ e~^* + o(l) , (3.12) 

where we use e~'^^* ' = e'^^ + o(l), which follows from 7"^ ^ T. The result then follows from 
Lemma 2.8 and Lemma 3.2. Items (iii) and (iv) are consequences of (i) and (ii) and the fact 
that 7r(S'°-) = 1/2 + 0(1). D 



(3.13) 
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From the previously stated results, one may conclude that there exists 5l, a decreasing se- 
quence tending to zero when L tends to infinity, such that for every t > 0: 

P^ [r°'- >t] ^e-^"*+<5L, 

P^^ [r°'- >t] ^e~^"*-(5L, 

TT{n+) ^ 1/2 - 6l 

7:+{S^'+)^1-6l. 

Given such a 5l, one defines S^'~^ to be 

^1,+ ._ ^0,+ n {ry e f^+ : P^' [r°'~ > t] ^ e"*^" - 3((5l)1/^ , Vt > } 
and 51'- = -S'^'+. 
Lemma 3.4. T/ie set S^'~^ satisfies 

^(51.+) = l+o(l), and 
^ ^2 ^ ^' (3.14) 

P'' [r°'" > *] = e"*'^° + 0(1) , uniformly inf^O, and r] £ 5*^'+. 

Proof. The lower bound in the second point follows from the definition of S"^' . For the upper 
bound, it is just a consequence of the fact that 

P« [r°'- >t] ^ P^ [r°'- > t] ^ e-^o* + 5l , (3.15) 

for any ^ G ^2, by monotonicity, where the last bound follows from point (Hi) of Lemma 3.3. We 
turn to a proof of the first point. For t ^ 0, one defines 

From the second line in (3.13), and using (3.15) for all £, G S"^'^'*: 

e"^«*-5L ^ P"^ [r°'" >t] ^ 7r+(5i'+'*)(e-T"* + 5L) + (l-^+(5i'+'*))(e-^o*-24/^). (3.16) 
This gives -26l ^ 2(5^/\+(5i'+'*)) - 261^^, i.e. 

7T+{S^'+'')^l-6l^\ (3.17) 

Next, define Si := S '~^' l 7o , i g N. We claim that for all L sufficiently large: 

S^'+ D f] SiH S^'+ . (3.18) 

1=1 

Indeed, let rj G ni=i '^i- ^i t ^ 6^ ^q /2, then clearly rj G S". r-i/4-| if L is large enough. 
In particular, 

ri[T^'- > t]^ P^'[rO'- > \tjoSl'^^^]Sl^^%^] ^ e--i^*-^T _ 2^^^ ^ e~^"* - 3(5y^ (3.19) 
If on the other hand t ^ S^ 7q" /2, then 

g-t7o_3ji/4^0^P''[r°'- >t], 
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provided 5l is small enough, i.e. L is large enough. This proves (3.18). Moreover one has, from 
(3.17) and the fourth line of (3.13), 

From the third line of (3.13), one gets that 

7r( nJS'''J S, n 5°'+) > (1 - 4/' - 6l){1/2 - 5l) • 

This last estimate together with (3.18) implies the first statement of the lemma. D 

The previous results allow us to compute the value of 71 = 751,-. 

Lemma 3.5. Let —71 be the largest eigenvalue of C ' , the generator of the process killed when 
it hits S"^'^ . We have 

7ir,ei = l/2 + o(l). 

It is important to recall that, in contrast to 71, the eigenvalue 70 of the process killed in C^ ' 
was estimated only up to a factor 2 (cf. Lemma 3.1). 

Proof. The inequality 7iTrci ^ 1/2 + o(l) comes from Lemma 3.4 and Lemma 3.1. Recall the 
definitions 

r-'i = inf{t^0,7?(t) g5^'-}, 

r+'° = inf{t^0,?7(t)G5°'+} . 

According to Theorem 1.5, (3.5), Lemma 3.1, Lemma 3.3(i) and Lemma 3.4, one can find a new 
sequence 5l going to zero such that 

zvf ((5°'+)^) ^ (1 - e-*/'^-')/2 + h for every t > , 

P5 [r+'° > t] ^ e-^«* - 5l for every t ^ and ^ G S^' , 

F^ [r~'^ > t] ^ e"^i* + 6l for every t^O. 

Next, define 

f+'° = inf {t ^ r-'i : r?(t) G 5°'+} 

to be the first time the process enters S^'~^ after entering for the first time S^'^ . One has 



2 



+ 6l^ i^i^iiS'^'+Y) ^ P^ [r^'i < t, r+'O > t] 



1{t-'^ < t}Fi^^ '') [t+'° > t]l ^ P^ [r-'i < t] min P*? [r+'° > t] 



^(l_e-7it_<5^)(e-t7o_^^). (3.20) 

We use this inequality for t = "/^ 6^ and get (using the fact that x — x^ ^ 1 — exp(— x) ^ x 
for X small enough) 

5l/V(27iTrei) > iSl^' - <5r - <5i)(l - 5^/^0/71 - Sl) - Sl. 
As 1 ^ 70/71 ^ 3 (for L large enough, cf. Lemma 3.1), all of this gives us 

7l^rcl ^ 



2[(1 - 61^' - ,5^/')(l - 351/3 _ 5l) - 61^'] ' 
which ends the proof. D 
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Once again assume that 5l is a sequence going to zero, this time such that one has for every 
t >0 

P^ [r^'- >t] ^e-^i* + (5L, 

P^^ [t^'~ >t] ^e-^i*-(5L, (3.21) 

7t{S^'+) ^ 1/2 - 5l. 

Note that the sets S^' are not yet good candidates for the sets S of Theorem 1.3, the reason 
being that (3.21) and Lemma 3.5 say that the hitting time of S^~ is exponential with the correct 
rate, but only if one starts from either the maximal configuration or from tt^ , while we want 
this to hold uniformly in the initial condition in S~^ . We need therefore a final step in order to 
fix this problem. We set 

52.+ := ^f] G 5!'+ : P^ [t^'~ >t] ^ e-^i* - 35^/^} , 

and define 5"^'" = — S'^'"'". The same computations of Lemma 3.4 prove 

Lemma 3.6. tt{S^'+) = 1/2 + o(l). 

Now we are ready to finish the proof of Theorem 1.3, with S^ := 5*^'^. Let —72 be the largest 
eigenvalue of the generator C^ ' of the process killed when it reaches S'^'". From Lemma 3.1, 
one has tt{S'^'~^) gap ^ 72 ^ 71- Therefore, Lemma 3.6 yields 

72 = (1/2 + 0(1)) gap. (3.22) 

Let T^'" be the hitting time of 5*^'". For any t] G S'^'~^ (this is actually true for any ry in 0), we 
get from monotonicity. Lemma 3.3(i), Lemma 3.6 and (3.22) 

r'{T^'~ > t) ^ P^(r2'- > t) ^ e-^2t + 0(1) + (2 - Att{S^'-)) ^ e~^i + o(l), 

where o(l) is uniform in t. On the other side, the definition of 5*^'^ and the obvious bound 
r^'" ^ r^'~ give that for any rj G S'^''^ 



>t)^ P''(t^'" > t) ^ e^^i* + 0(1) ^ e 2T,ei + 0(1), 



where the last inequality comes from Lemma 3.5. D 

3.3. Proof of Corollary 1.7. We use the same notation that in the previous proof. We set 
r := (S'2'+ U S'^'~y, and denote the local time spent by the Markov chain {r]s)s>o in T by 



Ht = Ht{T) := / l^(5)grds. 
Jo 



(3.23) 



Notice that, if vf^'^ := 7r(- | 6*2 ), then Hvr^'^ — 7r^|| = o(l). Equation (1.10) implies that there 
exists a sequence 6l going to zero such that 

— 1/2 

Integrating the first equation between zero and 6j^ Tj-ci one gets 



E*^'' 






< ^^/'Trel . 



We set 
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Using Markov's inequality we obtain 

Moreover, by Theorem 1.3, there exists a sequence S'^ going to zero such that 



-2,- ^ rn a1/16^ 1 , , re-l/2 



^2,- 



r^'- G [0, <5i/'Vrci] U [6l'^'Z,u oo)J < 5'l, Vr/ G 5' 

On the event {r^'~ G [SJ T^ei U 5 r Tlei]}, one has H^i,- < H i/2rp and hence, for every 
ri G 53'+, 

Tj \ xl/16 2,-1 / rl/S , X' 

D 

4. Mixing time upper bound 

In this section we prove Theorem 1.1. Our approach wiU also yield a proof of Theorem 1.2. 
The main ideas of the proof can be sketched as follows. 

Step 1. Lemma 2.8 shows that after a burn-in time 0{Lp'^'^) the distribution f^ has a smooth 
density w.r.t. the equilibrium vr. The first step consists in using this fact together with (2.6) 
and Lemma 2.4 to reduce the mixing time upper bound to a lower bound on the spectral gap 
of the chain. 

Step 2. To bound the spectral gap we decompose the polymer configurations using the variables 
(Tx = sigT].{rjx) introduced in (1.7). From the decomposition estimates in Proposition 2.1, we 
shall roughly obtain that the spectral gap of the chain is bounded below by gap^ x gap^ , where 
gap^ denotes the spectral gap of the heat bath dynamics for the variables a, while gap^ stands 
for the spectral gap of the polymer with a wall (i.e. the polymer constrained to be non negative). 
From [7], we know that gap_|_ = Q{L~'^). 

Step 3. To prove a lower bound on gap^ we shall perform a second decomposition, this time by 
fixing the number of crossings (i.e. the number of sign switches) in the configuration a. Another 
application of the bound from Proposition 2.1 will then show that gap^ is roughly bounded 
below by a product of two spectral gaps, say gap* and gap* . Here gap* = min„gap*„, 
where gap* „ is the spectral gap of the dynamics on the variables a constrained to have n 

(2) 

crossings, while gap* denotes the spectral gap of a birth and death chain associated to the 
number of crossings. We establish a lower bound gap* = i7(l). Moreover, we show that when 
n = 1, one has gap* ^ = i7(L~^'^). To prove a similar bound for every n we introduce a new 
dynamics involving a fixed number n of crossings: with rate 1, independently, each crossing 
equilibrates its position between the two neighboring crossings positions. If gap" denotes the 
spectral gap of this process, a comparison argument shows that gap* „ ^ gap* i x gap" . 

Step 4- The final step consists in obtaining the lower bound gap" = ri(L^^). The first ob- 
servation is that if n ^ elogL this estimate can be obtained by means of a direct coupling 
argument. The proof of the estimate for larger values of n is based on a block dynamics argu- 
ment which allows us to reduce the problem to the case of n ^ e log L crossings. The analysis of 
the block dynamics uses a further coupling argument. It is worth observing that the coupling 
arguments used here make crucial use of the heavy tailed nature of the distribution of excursions 
at equilibrium; see Lemma 2.2. 
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Before starting the actual proof, let us pause for a few remarks. The lower bound on gap^ 
described in Step 3 and Step 4 above is sharp (up to 0{L^) corrections). As detailed in Section 4.5 
below, Step 3 and Step 4 will essentially prove Theorem 1.2. On the other hand, the final bound 
Tmix = 0(L^'^^^^^) for Theorem 1.1 is likely to be off by a factor 0{L?). As explained in 
Step 1 above, this comes from the use of a decomposition estimate that involves the product 
gap^, X gap_|_ rather than the minimum min{gap^,,gap+}, as it would be the case if one could 
efficiently decouple the mode associated to the variables a from the rest. 

The following four subsections will develop the four steps described above in the given order. 
However, we warn the reader that, because of various technical obstacles, the above plan will 
not be followed very strictly and several detours will be needed. 

4.1. Reduction to spectral gap. We start with the implementation of Step 1. For later 
purposes it is necessary to consider a variant of the original dynamics which avoids (very unlikely) 
configurations with too many crossings or too many zeros between consecutive crossings. 
Call X the number of crossings in a configuration ry: 



x(^) 



X- 



L-2 

E 

=-L+2 



^iVx = 0, rj^^i / ■>], 



x+l) 



Note that only sites x G El \ ({— L} U {L}) appear in the summation. Define 7^ 

^L-2 

JX = 



(4.1) 

Hvx = 



0, Vx-i / Vx+i), so that X 



T.x=-L+2^x, and write 
^0 = --^ , ^x 



L. 



(4.2) 

and if 1 ^ j ^ X) l^t S,j denote the position in {—L + 2, . . . ,L — 2} of the j'-th "1" in the 
sequence {7_i_|_2, . . . ,7l-2}- Thus, Ci)---;?x denote the positions of the internal crossings. 
Finally, denote by N{^i,^i+i) the number of zeros in the path 77 strictly between ^j and ^j+i. 
See Figure 2. 




Figure 2. A configuration rj of the polymer with xiv) — ^ internal crossings in 
positions ^1, ^2, C3,'^4- Note that in this case Af(Co,Ci) = N{^2,i'i) = iV(C3,C4) = 
0, and A^(Ci, ^2) = ^{^A-, £,b) = 1- Below, the corresponding configuration of signs 
a = sign(?7). 



Fix a constant Cq > and define the event 

9.° = Irj eVt : xiv) ^ Co log L , and max N{^i, ^i+i) ^ Cq log L i . 

L 4=0,. ..,x J 



(4.3) 



Clearly rj G {0,°y implies that there are at least Co log L zeros in the path rj. Since A < 1, the 
bound of Lemma 2.2 shows that the number of zeros is exponentially integrable at equilibrium. 
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Therefore, for any p > 0, taking Co{p) large enough, we have 

TT{n") = 1 + O(L'P) . (4.4) 

The reason for introducing the restricted set ri° will be apparent in the sequel. For the moment, 
we point out that the restriction x ^ Cq log L is essential for our estimates in Section 4.4, while 
both restrictions x ^ Co log L and N{S,i,$,i-\-i) ^ CologL will be needed in the estimate of the 
parameter 7 appearing in the decomposition of the spectral gap; see the proof of Proposition 4.2 
and Proposition 4.4 below. 

Next, consider the polymer process restricted to stay in the set fi", i.e. the continuous time 
Markov chain with state space il", and generator (1.1), where the rates rx,±{rj) are replaced by 

Let also n" denote its reversible invariant measure, which is easily seen to coincide with 7r(- 1 17°). 
Let Tj.°[ denote the relaxation time of the process defined above. 

Lemma 4.1. Assume that T°^^ = 0{U') for some p > 0. There exists c > 0, such that for any 
e > and for all L ^ LQ{e): 

T,nix^^'+' + cr°ilogL. (4.6) 

Proof. Let t = L?"^^''^ + s, with s = ciT"^^ fo^ some ci to be fixed below. We prove that 

lkf-^ni^l/2. (4.7) 

From Lemma 2.4 and Lemma 2.5, (4.7) implies that Tmix ^ C2(-L^^^'^ + s) logL, for some other 
constant C2, which implies the lemma. To prove (4.7), we first introduce some notation. We 
write Vu" for the distribution at time u of the state of the Markov chain restricted to il°, when 
the initial configuration is some 7? G i7°. If 77 ^ il° we define Vu° = (5^ for all u. Next, we write 
Vu° = Sn/ti(^)'^«'° for a probability measure /i on Q.. Using symmetry we can write 

\\v'^ - z^/ll ^ 2||z^f - ufW + 2\\uf - uf'°\\ + \\uf^° - vf^°\\ . (4.8) 

We start with the observation that 

H - ^f II = \Ws^-' - <^ll ^ \Wts - vr+ll = o(l) , 



where the first bound is obtained by writing j^s'""' —i^g = Jii's — i^s )p(^) v') with p the maximal 
coupling of VfLg and vr'^, and the last bound follows from Lemma 2.8. 
We turn to the last term in (4.8). Let tt~^'° = tt^{- \ Q°) and observe that 

||7r+ - 7r+'°|| =0(1). 

This last bound follows easily from (4.4) and (2.27). Moreover, the bound (2.6) applied to the 
process restricted to 0° yields 

for some c > and for all u > 0. Therefore, the third term in (4.8) can be made smaller than, 
say, 1/4 by taking ci large enough in the definition of the time s. It remains to prove that the 
second term in (4.8) is o(l). Since the initial condition is sampled from the same distribution 

vr"'" we can couple the two processes {v^ )« > and (z/^ '°)„ ^ in such a way that they coincide 
until the first time when the unrestricted process exits from the set Q°. (Note that this time 
can be zero.) Thus, 

\Wf - ^fn\ ^ E ^+(%)lP(3n ^ s : z;2" ^ ^°) • 
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From (2.27) we know that vr"'"(?7o) ^ (2 + o(l))7r(r/o), so that the time-invariance with a union 
bound imphes 

\Wf - '^f'l ^ 2.L(1 - ^(0°))(2 + 0(1)) , 
where we use the fact that the average number of updates up to time s is bounded by 2sL. Since, 
by assumption, s = 0{U') for somep, we can use (4.4) to conclude that ||z^J —v^ '°|| = o(l). □ 

The fohowing three subsections wiU focus on the upper bound T"^^ — 0{L^''^^'^^'^). Once this 
bound is estabhshed, Theorem 1.1 wih fohow immediately from Lemma 4.1. 



4.2. Decomposing along crossings configurations. Recall the definition of the variables 
a G {—1, +1}*^^ given by ay = sign(77j^), where Ol is the set of sites in {—L, . . . ,L} with the 
same parity as L + 1. Note that the field a specifies uniquely the field ^ defined after (4.2), 
while ^ specifies o" up to a global sign switch; see Figure 2. The space O" can be decomposed 
into disjoint subspaces 

^° = ^aK , (4.9) 

where 0° denotes the set of ry G $7° such that sign(?7^) = ax for all x G Ol. Let So denote the set 
ofaUcr G {-l,+l}<^i such that il° / 0, i.e. theset of a G {-1,+1}'^^ such that x(cr) ^ Co log L, 
where x{')i defined in (4.1), is seen as a function of a = sign(ry). Consider the continuous time 
Markov chain on So with infinitesimal generator 

G^{a) = Y, U^) bK) - V'(^)] (4-10) 

x£Ol 

where 99 : 5o i— )■ M, cj^ is the configuration a fiipped at x, i.e. it is defined as a everywhere except 
at X where it equals —ax, and the rates 6x{(y) are given by 

dxia) =Y.7:{^\ K) [<,+ (^) l(a(r?^'+) = CJ-) + rlAv) M'^iv^n = cx^)] • (4.11) 

The rates r° _|_ are given in (4.5). Note that the measure 

,^oi^)= Y,<^\^°), (4-12) 

is the reversible distribution, i.e. i'o{cr)6x{a) = fo(o"^)^a;(cr^) holds for all x G Ol and a G So- 

In words, the process with generator Q is described as follows. Attach independent rate 1 
Poisson clocks to all sites x G Ol. Let a be the current configuration. When site x rings, choose 
a configuration 77 sampled from the distribution 7r(- 1 Q°) and set r]' = 77^'"*" with probability 
rx^if]), rj' = ri^'~ with probability r^_{r]), and rj' = rj with probability 1 — r° _^_{r]) — r°_{r]). 
Finally, update a to a' given by a'^ = sign(r/^). Let gap"^" denote the spectral gap of this Markov 
chain. 

Proposition 4.2. There exists c > such that for all L: 

T°i^cL2(logL)3(gap'5o)-i. 

Proof. We apply Proposition 2.1 with the decomposition (4.9). To each a we can associate the 
continuous time Markov chain with state space i7°, defined by the generator (1.1) with the rates 
rx+irj) replaced by 

<±(r,)=<±(r,)l(r/-'± G K) , (4-13) 

where the rates r° ±{ri) are defined in (4.5), with reversible equilibrium measure tt°{- \ r2°). Call 
gap*^ the spectral gap of this Markov chain. For a given o", this corresponds to independent 
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continuous time Markov chains for each interval {^i, . . . ,.^j+i}, where the crossing positions ^i 
have been defined in (4.2). On a given interval {S,i, ■ ■ ■ ,^i-\-i}, we have a polymer dynamics 
with a horizontal wall constraint (polymer above or below the wall depending on the sign of the 
field a inside that interval). Moreover, within each interval the polymer is constrained to have 
smaller than Co log L zeros. Let gap'^'* denote the spectral gap of this process on the interval 
{^i, . . . , ^i+i}. From the independence recalled above, one has 

gap'^ = min gap'^'* , (4-14) 

i=0,...,x 

where x is the number of interior crossings defined in (4.1). It follows from Lemma 4.3 below 
that (gap*^)"^ = 0(L^ log L), uniformly in o" G So- 

Next, observe that the generator Q in (4.10) coincides with the generator C from Proposi- 
tion 2.1 for the present choice of the decomposition. We can then use the bound of Proposition 2.1 
with A = gap-^", and Amin = min^g^^ gap'^: 

r°i ^ c (gap-^-)-! max{l, L^ logL(gap'5° + 7)} , (4.15) 

where c > is a constant, and 

L-l 



. maxmax 2^ [r°_+(r,) l(a(7?-'+) = a^) + r°,_(r/) l(a(ry--) = a^)] . (4.16) 

It is immediate to check that gap ° is smaller than a constant, so that T°^^ ^ cL^ log L(l + 
7) (gap °)~^, by (4.15). It remains to give an upper bound 7 = 0((logL)^). From the definition 
(4.16), we see that 7 is bounded above by the maximum over r] € Q° of the number of sites 
X £ Ol such that r]x~i = Vx+i = 0. By definition of the set i}°, the latter quantity is bounded 
by (co log L)^ in our setting. This ends the proof. D 

We turn to the lower bound on the gaps defined in (4.14), that was needed in the proof of 
Proposition 4.2. The bound (gap"")^^ = 0{L^ logL), uniform over a G So, is an immediate con- 
sequence of Lemma 4.3 below and (2.8). Consider the polymer dynamics under the constraints 
N{r]) ^ M (where N{r]) is the total number of zeros and M is a positive constant) and 77 ^ 0, 
i.e. let r_|_^jv/ denote the set 

r+^M = {?? G rj : N{ri) ^ M , and r/^, ^ for all x} , 

and write p^ for the law at time t of the corresponding process. This is the continuous time 
Markov chain with generator (1.1) with rates replaced by rx,±{ri)l{r]^' G T-^-^m) and with 
reversible measure 7r(- | T^^m)- 

Lemma 4.3. There exists a constant c > such that for any rj G T^^m and for any M ^ c logL.' 

IIp?^ - 7r(- 1 r+,M)|| = 0(1) , T = cLHogL. 

Proof. Let ^^ denote the evolution without the constraint N{r]) ^ M, i.e. the Markov chain 
with generator (1.1) with rates given by rx^±{ri)'^{ri^'^ ^ 0) and with reversible measure tt"' = 
7r(- 1 ry ^ 0). The mixing time of this "wall" constrained model has been analyzed in [7, Theo- 
rem 3.1], where it is shown that, for some constant c > 0, for all initial r/ ^ 0: 

||^?-vr-KcL2exp(--^). (4.17) 

As in (4.4), standard equilibrium estimates imply that for any p > {) one can choose c > such 
that for all M ^ c logL: 

7r»(iV >M) = IK"- - 7r(- I r+,M)|| = 0{L-P) . (4.18) 
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Next, observe that by monotonicity fil is stochastically dominated by p^. Let t"^ denote the 
hitting time of the set A''(-) > M for the process with law (/U^)t 5. 0, and introduce the event 

Note that, for fixed t, Gt is an increasing event, and therefore p^l{Gt) ^ pl{Gt), for any s, t ^ 0. 
Let us fix now t = ci L'^ logL. For any p > 0, using the time-invariance and a union bound 
with the fact that the expected number of updates up to time t is bounded by 2tL, one has 

^7r"'(r?)P(r'' ^ t) ^ 2iL^"'(iV > M) = 0{L-p) , 

V 

where we use (4.18). Markov's inequality then implies that 

TT'^iGf) = 0{L^-P) . 
From (4.17) we then deduce that if ci is sufficiently large (in the definition of t): 

p^Gt) ^ nUGt) = 0{L'~n . (4.19) 

On the other hand, writing p2ti^) ~ 7r'"(yl) = ^ , p1{ri'){pf (A) — ^^{A)) for any set A one has 

\\pi-^n\^ ^ pKv'Wt - ^n\ + pnGt) . 

v'eGt 

To estimate the first term above, note that the processes with laws {p^)t^o,{pf)t^o can be 
coupled in such a way that they coincide until time t'^. Therefore, by definition of Gt, for any 
V' G Gt: 

\\p1' -tt-^W^ ¥{t'>' ^t) + \\p1' - vr^'ll = 0{L-^) , (4.20) 

where we have used again (4.17) to bound the last term above. In conclusion, using 

WpI - 7r(. I r+,M)|| ^ llvr- - 7r(. | r+,M)|| + \\pI - tt^W , 

together with (4.18), (4.19), and (4.20), we arrive at \\p2t — vr(- | F+^m)!! = 0{L^^) which implies 
the desired estimate. D 

Thanks to Lemma 4.1 and Proposition 4.2, Theorem 1.1 will follow from the estimate 

Proposition 4.4. For any e > and for all L ^ -^o(^)-' 

gap-^" ^ L^i"^ (4.21) 

The following two subsections are devoted to the proof of Proposition 4.4. 

4.3. Decomposing according to the number of crossings. We first decompose So accord- 
ing to whether the first excursion has positive or negative sign, i.e. 

5o = 5+ U 5- , S^ = {a£So : a_L+i = ±1} • (4.22) 

An application of Proposition 2.1 with the decomposition (4.22) yields 

gap-^" ^ min{A/3, ~XXmin/{X + 87)} , (4.23) 

where Amin = gap"^ denotes the gap of the process restricted to S~^ (by symmetry, this equals 
the gap of the process restricted to S~), while A is the gap of the symmetric two state Markov 
chain with transition rate 

c(+,-) = c(-,+) = V i^o(cr|5+)6l_L+i(cr), 7= max6'_L+i(o-), 
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9^L^i{a) being the probability of a flip at x = —L + 1 in c; see (4.11). The probability measure 
Uo is defined in (4.12). Clearly, 7^1. Since at equilibrium the polymer has a uniformly positive 
probability of taking the value ??-l+2 = 0, cf. (2.24), it follows easily that c(+,— ) ^ c(A) > 
for some L-independent constant c(A). Thus, A is of order 1, and for some c = c(A) > one 
obtains 

gap"^" ^ cgap-^^ . (4.24) 

Next, consider the number of crossings x defined in (4.1). Since x is a function of the signs a 
only, we can write xi^^) ^oi the number of crossings in a given a £ So- Thus, the space S~^ is 
partitioned as 

5+ = U-=o'5+'", (4.25) 

where m = [co log LJ, and 5+'" = {a £ S^ : xi^^) = "-}• We apply Proposition 2.1, this 
time with the decomposition (4.25). Thus, gap can be bounded below as in (4.23) where 
Amin = mino <in !^m gap"*"'" now stands for the minimum of gaps of the process restricted to 
5+'", while A stands for the gap of the random walk on {0, . . . , m} with transition rates 



c(n,n±l)= Yl E '^°(^l'5^'")^-(^)l(^"^'5^'"^')' 


(4.26) 


cr(^S + '"xGOL 




c{n,n±2)= Yl Y^o{a\S+ne^{a)l{a-£S+'^^'). 


(4.27) 


aeS + '-^x&OL 





Note that the transition n — )• n + 1 is necessarily obtained by creating a new crossing at the 
rightmost position ^„ = L — 2, while the transition n — )■ n + 2 can be obtained by creating two 
adjacent crossings between two existing consecutive crossings. These are the only transitions 
that are induced by the single flips o" — )■ cj^. By construction, the reversible invariant measure 
of the chain defined by the rates (4.26), (4.27) is given by 

n{n) := Vo{x = "^ I S^) = ^^ix = "^ I ^°) : n = 0, . . . , m . 

Lemma 4.5. The gap A of the chain defined by (4.26), (4.27) satisfies c ^ A ^ c^^ for some 
c = c(A) > uniformly in L ^ 2. 

Proof. By Lemma 2.2, tt(x = it-) is exponentially decaying in n, and the same applies to ^(n) = 
^(x = ^ I ^°)) up to 0{L^P) corrections; see (4.4). To prove an upper bound on A one can take 
the test function x counting the number of crossings in the variational principle defining A. The 
variance of x w.r.t. fi{-) is a positive constant. On the other hand the Dirichlet form can be 
bounded from above by 

m~l 

4 Y /"H[c(f^, n + 1) + c{n, n + 2)] . (4.28) 

n=0 

Observe that for every n, and every a G 5^^'", one has 

xeOL 

Indeed, each excursion can only contribute 0(1) to this sum since at equilibrium, in the delocal- 
ized phase, the expected number of zeros between consecutive crossings is finite (depending on 
A < 1); see Lemma 2.2. Now, (4.29) implies that (4.28) is bounded above by c ^^ /u(n) n = 0(1). 
This ends the proof of A ^ c for some constant c. 

To prove a lower bound on A we can neglect the additional rates (4.27). Let us check that 
the rate c(n, n + 1) in (4.26) is bounded away from zero uniformly inO^n^m— 1, L^2 {m 
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being the maximal number of crossings allowed) . We have 

c{n, n + 1) ^ c i/o(^„ = L - 4 I 5+'") , (4.30) 

for some constant c = c(A) > 0, where ^„ denotes the rightmost internal crossing. We introduce 
some extra notation to characterize more explicitly the measures Vn^o '■= ^o{' \ S~^'^). Let 0.- '° 
denote the set of polymers 

^T° = {^ ^ • Vo = Vj = 0, rjx+i = r]x±l, x = 0,...,j -1, and N{t]) ^ Co log L} . 

Note that these paths start at and end at j, so that j must be even for this set to be non 
empty. Define the probability 

where z+^o is the normalization. From Lemma 2.2, (2.23) one has that 2^^ Z- '° = 0(j^^") 
for large j, so that the probability p+^o is well defined. Then, it is not hard to check that the 
measure Un,o introduced above is given by 

p+,oixi + L)p+^o{x2- Xl)--- p+^Xn- Xn~l)p+,o{L- Xn) /,oi\ 

Un,o{xi,...,Xn) = ^ (4.31) 

where zy„,o(xi, . . . ,x„) = f„,o(6 = xi,...,S,n = x„), if -L < xi < • • • < x„ < L is any 
allowed configuration of the crossing positions ^i, . . . , ^„, and Z° ^ is the normalizing constant. 

In particular, Un^o is a product measure p_|_ " over n + 1 positive increments Xj+i — Xi (with 

xo = —L,Xn+i = L), conditioned to have X]"=o(^*+i ~ ^«) ~ '^^■ 
Going back to (4.30), with the notation in (4.31) we have 

p+,o(4)/.5r^(Er=ite-^.~i) = 2L-4) 



l^oiCn = L-4\S^ 



pa"^'^(Er=^i'(6-6-i) = 2L) 



Since vL ^ m ^ n, and p-|-^o(4) > uniformly, with the same arguments of Lemma 4.12 below 
one easily sees that I'oiCn = L — 4|5"'"'"') > uniformly. (We omit the details here to avoid 
repetitions) . 

Once we have that the rate c{n^ n + 1) in (4.26) is bounded away from zero uniformly, a 
standard bound using the isoperimetric constant and Cheeger's inequality (cf. e.g. [12]) implies 
that A is bounded away from zero uniformly. Alternatively, the same conclusion can be derived 
from more refined bounds for birth and death chains [17]. D 

We return to the application of Proposition 2.1 with decomposition (4.25). The constant 7 is 
now given by 

7 = max max V 6 Jo) [l((J^ G 5+'"^^) + 1(ct^ G 5+'"^^)! _ 
xeOL 

Recall (4.29). Since n ^ m = O(logL) we obtain that 7 = O(log-L). In conclusion, using (4.24) 
and Lemma 4.5, we arrive at 

gap"^" ^ c(logL)"^ min gap+'" . (4.32) 

1 ^ n ^ m 

Thanks to (4.32), the proof of Proposition 4.4 has been reduced to the proof of a lower bound 
gap^'" = ri(L^^'^^^), uniformly in n ^ Co log L. Note that the case n = does not appear 
in (4.32). Indeed, in that case the corresponding equilibrium is concentrated on the single "all 
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plus" element 5"*"'". We start with the simplest case, i.e. n = 1, which can be analyzed by rather 
standard arguments. 

Lemma 4.6. There exists a constant c = c(A) > 0, such that 

gap+'^ ^cL^^/^ (4.33) 

Proof. When n = 1, the process restricted to 5^'^ consists of a nearest neighbor random walk 
on El, the sites in {—L + 2, . . . , L — 2} with the same parity as L. Letting cr{x) denote the 
element of 5"*"'^ with unique crossing at x, the corresponding jump rates are given by 

c(x, x±2)= Y^ ey{(7{x))l{a{x)y = a{x ± 2)) . (4.34) 

These rates are of order 1 (except at a; = — L + 2 where c(— L + 2,— L) = and at a; = L — 2 
where c(L — 2, L) = 0) since there is a uniformly positive probability for the polymer conditioned 
to have signs cr{x) of having a at a; ± 2. Moreover, the reversible invariant probability measure 
p{x) for this chain is given by 

p{x) = 7r"{r,^^i,r,^, r/^+i = (1, 0, -1) | sign(r/) G 5+'^) . (4.35) 

Recalling (4.4) and Lemma 2.2, the event sign(r/) G S^'^ has uniformly positive probability at 
equilibrium, and p{x) x po{x), where pQ is the probability 

Po{x) = ^-^^ ^ , x€El, (4.36) 

where zl is the normalizing constant. Since the rates (4.34) are of order 1, a standard comparison 
argument shows that it is sufficient to prove the bound gapg ' ^ c L~^ " for the gap of the 
chain with reversible measure po given by (4.36) with rates defined by co(x,x + 2) = 1 and 
co(x, X — 2) = po{x — 2)/po{x). The latter process has the Dirichlet form 

L-4 

£o{f,f)= Yl Po{x)[f{x + 2)-f{x)]\ (4.37) 

x=-L+2 

On the other hand, writing f{y) — f(x) = Yl^j=xifO + 2) — fU)], x < y, and using Schwarz' 
inequality, the variance Var^o (/) can be bounded above by: 



^poU ) — 2 Z^'"" 



Varpo(/) = - > ^Po{x)poiy)[f{y) - fix)Y 

po{x)po{y) 



x,y 

L-4 



From (4.36) one has 



^2L Y Po(j)[/(j + 2)-/(j)PEE 

j=-L+2 x^j y>j 



^,^YypjMp^ = o{lI: 



Poij) 



' .^.fe p'^'^ 



It follows that gapQ ' ^ cL ^'^, which implies (4.33). D 

Remark 4.7. The bound of Lemma 4.6 is optimal in the sense that 

gap+'i ^ c-i L-^/2 _ (438) 

This can be seen by taking a test function ip{x) in the variational principle for gap^'^, of the 
form f{x) = g{x/L), where g : [—1,1] i— )■ M is given by g{s) = — 1 if s < —1/2, g{s) = 1 if 
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s > 1/2, and is linear between —1/2 and 1/2. With this choice one has that the variance of (p 
w.r.t. p defined by (4.35) is of order 1, while the Dirichlet form - given by (4.37) with p instead of 
po - is bounded above by the probability of x G [— L/2, L/2] times L~'^, i.e. L"^'^. This implies 

(4.38). 

We turn to the proof of a lower bound on gap^'" for n > 1. Recall the definition (4.31) 
of the measure i'n,o- We introduce a further dynamics. We view the ^i as particle positions. 
Each particle i = 1, . . . ,n has an independent Poisson clock with parameter 1. When particle i 
rings, we freeze all positions ^fc, k ^ i, and update ^j with the new position ^'- sampled from the 
conditional distribution f^ := i'n,oi: \^k, k ^ i). The Dirichlet form of this process is given by 

n 

Cq°(/, /) = E ^->o [Var,. (/)] , (4.39) 

where Varj^i (/) = ^'^(/^) — J^nif)'^ is the variance conditioned on the values of ^fc, k ^ i. Let 
gapc4° denote the associated spectral gap: 

gcT(/,/) 

TVar,„,„(/)' 



gap^° = irif "^;;V VV, , (4.40) 



where the infimum ranges over all functions of the crossing positions ,^i,...,^„. The next 
estimate allows one to reduce the proof of Proposition 4.4 to the proof of a lower bound gapeq° = 
VL[L^^), for every n ^ CglogL. 

Lemma 4.8. There exists a constant c = c(A) > 0, such that for all n^ 1. 

gap+'"^cL-5/2gapSi°- (4.41) 

Proof. For n = 1 this coincides with the result of Lemma 4.6, since in this case i^i^o[Var^,i(/)] 
coincides with Varj^^ „(/). The general case follows from the observation that Lemma 4.6 can be 
applied with the segment {^j_i, . . . , Ci+i} replacing the usual {—L, . . . , L}, to obtain 

Var,, (/) ^ - (e.+i - C^-lf/'£\f, /) ^ - L5/2r(/, /) , 

" C Ci 

where ci = c/2^'^, and £^{f,f) denotes the Dirichlet form of the nearest neighbor random 
walk on El, corresponding to the rates (4.34) (with the segment {S,i-i, ■ ■ ■ ,S,i+i} in place of 
{— L, ...,L}). Taking the z^n^o-expectation, and summing over i = l,...,n one obtains the 
estimate 



1 



t=l 

The desired conclusion follows from the observation that Y17=i ^n,o [^*(/) /)] = £^'^{fj /)) where 
£~^'^{f, f) is the Dirichlet form of the process restricted to S~^'^ with spectral gap given by 
gap+'". D 

4.4. Lower bound on gapcq°. To complete the proof we need the lower bound gapeq° = Vl{L^^). 
We first remove the restriction N{rj) ^ Co log L in the definition of the measure Vn,o- Namely, 
introduce the probability measure 

/ X _ P+{XI + L)p+{X2 - Xi) • • • p+{Xn - Xn~l)p+{L - Xn) , . 

l^n\Xl, . . . , Xn) — „ ) v^-^^j 
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defined by the kernel 

P+(J) = ^^ , Zi=y^ A^('') , (4.43) 






where 

^/ = {?] ^ : r/o = r/j = , jy^^+i = Ty^^ ± 1 , x = 0, . . . , j - 1} . 

Note that (4.42) coincides with (4.31) except for the removal of the constraint N{r]) ^ Co log L in 
the definition of the kernel p+. As in (4.4) one can check that p+ is equal to p+^o up to 0{L~p) 
corrections for some large p > and therefore i^nixi, ■ ■ ■ ,x„) = i^n,oixi, ■ ■ ■ ,a;„)(l + 0(L~p)), 
uniformly. Proceeding as in (4.39) and (4.40) we define 

n 

Cq(/,/) = E^"[Var..(/)] , (4.44) 

and gap" , the spectral gap associated to the measure f„ and the Dirichlet form (4.44). From 
the previous observations we see that Yar^^{f) ^ c Vart,^^(/) and Scq'ifif) ^ CiS" (/,/), for 
some constant c > 0, for every function /. Therefore gap" ^ c^ gapcq°. For later purposes it is 
important to keep track of the dependence on L in our notation, and therefore we write gap" (L) 
below. The desired bound gapcq = J7(L~^) follows from the comparison mentioned above and 

Proposition 4.9. For any e > 0, there exists LQ{e) > such that for L ^ Lq and n ^ Co log L 

gap^q(L) ^ L-^ . (4.45) 

As a preliminary step towards the proof of Proposition 4.9, we establish the following lemma. 

Lemma 4.10. There exists c > such that for any n, uniformly in L '^ n: 

gap:q(L)^ce-"/^ 

The crucial point of the above estimate is that it does not depend on the size of the system 
L, but only on the number of particles n. Note that Lemma 4.10 gives the required lower bound 
(4.45) if we take n ^ ei log L with ei suitably small. The case of larger n will be considered 
afterwards. 

Proof of Lemma 4.. 10. Every particle has an independent Poisson(l) clock. When particle i 
rings we update its position according to the equilibrium measure conditioned on the neighboring 
positions ^i-i and ^j+i, that is the probability for the i-th particle to be in position ^[ = ^i_i + /c, 
/c = 2, 4, . . . , ,f j+i - ^i_i - 2, is given by 



P+{k)p+{ii+i -ii^i -k) 



ii+l 



-6-1-2 






where Z= >^ p^(^k)p+{^,+i - ^,^1 - k) . (4.46) 



Z ' ^ 

k=2 

An important property of (4.46) is that the probability a of the event ^'^ = ^j_i + 2 (or ^'- = 
^i+i — 2) is positive, uniformly in the length ^j+i — ^j_i. The idea is to use this property to 
prove that with probability at least 2'^"'^ uniformly in the initial configuration, the process hits 
the minimal configuration 

6 = -L + 2, 6 = -i + 4, . . . , C„ = -L + 2n , (4.47) 

before time n?. Once this result is available, one concludes with a simple coupling argument. 
Indeed, using an independent coupling, the total variation distance between two evolutions at 
time t = mx n?, with arbitrary initial conditions, is bounded above by (1 — q^"/4)™ ^ g-ma /4_ 
Thus, the mixing time of this chain is at most 8n^a~^". Using e.g. (2.8) we obtain the lemma. 
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To prove the above claim, we use the notation {Tm,^m)m^ i for the sequence of updating 
marks, that is the m-th update occurs at time Tm and it concerns the £m-th particle, where 
im is a number in {!,••• ,n}. Consider the event E that there exist integers 1 ^ 0(1) < 
0(2) < • • • < (j){n) such that t^(„) ^ n^, i^i) = i, and such that for every i = 1, • • • , n and 
0(i — 1) < j ^ 0(«), one has tj / ^^(i-i) (where 0(0) = and £o = 0). In words, E is the 
event that within time -n? there has been a sequence of n updates at times T(^(i), • • • , T(^(n), such 
that the update at time TMi\ concerned the i-th particle and such that the i-th particle is not 
touched again before time T^u^iy Conditioned on the event E one has a probability of at least 
a" of hitting the configuration described in (4.47). Indeed, at time T(^(i) we set ^i = 2 with 
probability a, at time t^(2) we set ^2 = 4 with probability a, and so on. Therefore, to prove 
the lemma it remains to show that the event E has probability at least 1/2. This can be easily 
seen as follows. Consider the event F that a sequence of integers 1 ^ V'(l) < ^"(2) < • • • < 'ip{n) 
exists such that r^(n) ^ n^, and £^(j) = i, for every i = 1, • • • ,n. Let t^^' denote the first time 
particle 1 is updated and define recursively r'*' as the first time after r^*"^-* when particle i is 
updated. Clearly, F = {r^") ^ n^}. Using E[t''"'] = n, by Markov's inequality the probability 
of F is at least 1 — n/'n?' ^ 1/2 for n ^ 2. On the other hand E = F, since one can define the 
sequence from ip by choosing 0(n) = ip{n), and for 1 ^ i ^ n — 1, 

0(n — i) = max{j < ijj{n — i + 1) : lj=n — i}. 

This ends the proof of Lemma 4.10. D 

Proof of Proposition 4-9. The proof of Proposition 4.9 is based on a block-dynamics argument 
that allows one to reduce to n ^ eilogL particles, in which case the result will follow from 
Lemma 4.10. Fix an integer K ^ n. A block is a collection of particles with adjacent labels, and 
our particles will be partitioned into A := [n/K\ non-overlapping blocks, in such a way that 
the first A — 1 blocks contain K particles, and the last block contains at most K — 1 particles. 
For the sake of simplicity, since it does not change any of our estimates, we will suppose that 
all blocks have exactly K particles, i.e. n = KA. With this notation, the configuration of the 
i-th block can be described by the variables 

i^K{^-l)+l,■■■,CK^), i = I, . . . , A . (4.48) 

As usual, the S,i are interpreted as particle positions or crossing positions. To define the block- 
dynamics, we consider independent Poisson(l) clocks on each block, when one of them rings we 
put all crossings of the relative block simultaneously at equilibrium conditioned on the position 
of all crossings belonging to the other blocks. That is, if Bi denotes the i-th block (4.48), and 
Var^. is the variance with respect to 

then the Dirichlet form of the block-dynamics is given by 

A 
£m,f) = Y.^n[ys.TBAf)], (4.49) 

for any function / of the particle positions. Call gapjjj the associated spectral gap. The gap of 
the original dynamics for a single block Bi is given by gap^(£), where i = CiK+i — ^{i-i)K is the 
length of the portion of the system occupied by the K particles in the i-th block: 



gape^(^) VarB,(/) ^ ^n [Var,. (/) | Bj, j / 
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where Var^^i (/) is as in (4.44). Using this estimate in (4.49), taking the i^n-expectation and 
summing over i one obtains 

gap^q(L) ^ gap^i X mmgap^(^) . (4.50) 

Note that this is a special case of a well known estimate that controls the gap of the original 
dynamics in terms of the gap of the block-dynamics (see e.g. [16]). From Lemma 4.10 we know 
that gap^(^) ^ e~'^ uniformly in the length £. Therefore, Proposition 4.9 follows directly 
from Lemma 4.11 below by taking K = eilogL, K ^ n ^ Co log L, with ei sufficiently small 
(depending on e) and A = n/K (the case n ^ ei logL being contained in Lemma 4.10). 

Lemma 4.11. For all A ^ 1, there exists c(A) > such that for all L and A ^ n ^ Co log L, 

gapSi ^ c(A) K-*^ , K = n/A. (4.51) 

Before proving Lemma 4.11, we need to establish some preliminary facts. Consider u^^ '■= i^m 
the equilibrium measure (4.42) of the system of n particles over the segment {— L, ...,L}. 
Passing to the increment variables Q := ^i — £,i~i, i = l,...,n + 1 one writes v^ as the 
conditional probability p^ ( • | Y17=i Ci = 2-^)- 

Lemma 4.12. For all n = o{L?-'^): 

i.f(Ci^2L-LV3) = _l_(i + o(l)). 



Proof of Lemma 4.12. Observe that 

.f (Ci ^ 2L - L'l^) = ^^^2L-I^^^P^i^)f^(^^C. = 2L-^) 

Recall that p+{j) ~ C-^-j~^''^, for some known constant c+ > 0; see (2.23). In particular, 

p+{i) = p+{2L) (l + 0(L-2/3)) , i G {2L - lV3, . . . , 2L} . (4.52) 

Also, since 2L ^ (n + 1)^, using [8, Theorem A], one has 

n+l 

^0{n+i) (^ ^ ^_ 2l) = (n + 1)(1 + o(l)) p+{2L) . (4.53) 

Therefore, 

n 

(n + l).f (Ci >2L- L'/') = (1 + 0(1)) p?") ( J] e [0, lV3)) . 

It remains to show that />+ (X^"=i Ci £ [0,^"^'^)) = 1 + o(l). However, this is an immediate 
consequence of the assumption n = o{L^'^) and well known estimates for heavy tailed random 
variables. D 



Proof of Lemma 4. 11. As a corollary of Lemma 4.12, one has that for L large enough, uniformly 
in 21 ^ Li/3 and in n = o{i^/^), 

ul\Ci >2i- Li/3) ^ vl\c^ >U- i^l^) ^ ^^^^ . (4.54) 
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We shall use this observation together with a coupling argument. Start the dynamics at two 
arbitrary initial configurations ^,^' of n particles. We call S,{t),S,'{t) the corresponding states of 
the dynamics at time t. It will be shown that for a suitable coupling P: 

P(e(l) = ^'(1)) ^ c(A)(5(i^ + 1))-^^ . (4.55) 

Once this estimate is available the conclusion follows easily. Indeed, (4.55) implies that at time 
T the total variation distance between S,(T) and ^'{T), is bounded above by 



:i-c(A)(6(i^ + l)) 



~4A\T 



^ e 



-TciA)(5{K+l))- 



and therefore the mixing time of the chain is bounded by 2c{A)^^ {5{K + 1))^ . Using (2.8) 
we obtain the same bound for the inverse spectral gap. After adjusting the constant c(A), this 
proves Lemma 4.11. 




Figure 3. Illustration of the coupling used in the proof of Lemma 4.11 in the 
case A = 2, n = 6 and K = 3. At time ti, all particles in the second block are 
placed to the right of L — L^'^. At time t2, all particles in the first block are 
matched to the left of —L + L^'^. At time t^, all particles in the second block 
are matched. 



To prove (4.55) we proceed as follows. We have A independent Poisson(l) clocks at each block, 
and we note {tm,Cm)m > i the sequence of update marks: for any m '^ 1, tm is the time of the 
m-th update and c^ G {1, • • • , A} is the label of the block to be updated at time tm- Consider 
the event E that within time t = 1 a sequence of updates ti, . . . , tA-i, ^A, ^A+i, • • • , i2A-i has 
occurred, such that Cm = A — m + 1 for tti, = 1, . . . A, and ca+j = i + 1 for i = 1, . . . , A — 1. 
Clearly, E has a positive probability c(A) > 0. On the event E we define a coupling of the 
two evolutions by using the same marks {tm,Cm), i-e. at time t^ we update block Cm in both 
configurations. We refer to Figure 3 for a representation of the case A = 2. At time ti we have 
to update the particle positions in the block B/^. The corresponding equilibrium is of the form 
L — S,K(A-i)- Then, the new particle positions will satisfy 

^Ki A-i)+i ^ L — L ' , with probability at least — -— — — - . (4.56) 

2[K + 1) 



i/f , with 2i 
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Indeed, either ^x(A-i) ^ ^ — ^^ and the requirement is automatically satisfied because of 
^ii'(A-i)+i ^ C_ft:(A-i)) oi' '^x(A-i) < L — L}'^, in which case 2£ ^ L-^'^ and the claim follows from 
(4.54). Therefore, using an independent coupling, one has a probability at least {A{K + 1)^)~^ 
for the event 

min{^;^(A_i)+i(ii),^x(A-i)+i(*i)} ^ L - L^'^ . 
Then, one updates the block -Ba-i at time t2, and so on until one updates the block B2 at time 
tA-i- Iterating the argument given above, one has a probability of at least (4(i^+ 1)^)~ "''^ for 
the event 

min{ei^+i(tA-i), eWi(*A-i)} ^ i - (A - l)Li/3 . (4.57) 

Next, at time tA, one updates the block Bi. Let us show that on the event (4.57) there is 
a coupling P of the two equilibria on block Bi (conditioned on the value of ^K+i{tA-i) and 
^ii-+i(iA-i) respectively) such that the event 

Ml := te(tA) = ^-(^a) ^ -L + Li/3, i = l,...,K}, 

has probability 

¥{Mi)-^{b{K + lfy^ . (4.58) 

Clearly, it suffices to show that £,K{tA) = C'xi^^) ^ — L + L^'^ with at least that probability 
(all other particles S,i{tA), ■ ■ ■ ,S,K-i{tA) are then automatically matched using the diagonal 
coupling since Bi is the first block and Co = Co ~ ~^ ^^ fixed). Setting 2i = C-ft'+i(tA-i) + L 
and 2i' = C/4-+i(*A-i) + L, we need to couple the measures i^'^iCK = i), ^"k {^'k ~ ^)- Consider 
first the problem of coupling i'i{i) := J^^i^K = -^ + ^1^-^^ — L + L^>^) and 1^2(0 := I'x i^K ~ 
L + il^K^ -L + LV3). We have 

\\vi-V2\\ = 0{L-'^'^). (4.59) 

Indeed, first note that, by (4.57), one has 

2L^2£, 2i' ^ 2L - (A - l)L^/^ . (4.60) 

With the notation p_,_ {i) := p^ w2j=i Cj = ^) one has 

Li/3 
\\vi -V2\\ = 2 IZ l^i(0-'^2(i)l 

2 h ' E™=i pf ^(m)p+(2^ - m) E™=1 />?^(m)p+(2f - m) 

Using (4.52) and (4.60) one has p+(2f-i) = p+{2i-j){l+0{L-'^/^)) uniformly in ^ i,j ^ L^/^. 
Therefore (4.59) follows from 



Ll/3 Cg)(i^)/.N /o/C 



-2/3^1 y \p p+ vvp+v^^-^; ^ n(r-^/^\ 
fiESL{pt''im)p+i2i-m) 

We turn to the proof of (4.58). We define the coupling of i^'^iCK = i), i^k H'k ~ *) ^^ follows. 
Flip two independent coins with head probability p = i^j^i^K ^ —L+L^'^) andp' = uj^ (C^ ^ — 
L + L^'^). If both coins end up being head, then sample the pair ^k^Ck according to the 
maximal coupling of z^i, 1^2 defined by ||i^i — J^2||- If the p coin is head and the p' coin is tail then 
sample independently S,k according to f|^(- 1 Cftr ^ — L + L^'^) and Cx according to u]^ {'\Ck ^ 
—L + L^''^). Similarly, if the p coin is tail and the p' coin is head then sample independently 
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^^ according to ^'^'(- \^'k ^ — L + V-/^) and ^k according to v'j^{- \ ^k > -L + L^/^). Finally, 
if both coins are tail then sample independently ^k,^'k according to z^|?(- 1 ^k > —L + L^'^), 
^K (■ I ?i4' ^ ~^ "I" L^ )■ Using this coupling, we obtain that 

P(Mi) ^ pp' (1 -\\v^-v^\\)-^ {A{K + lf)-\l + 0(1)) , 

where we have used the fact that niin{p, p'} ^ 1/2{K + 1) by (4.54). This ends the proof of 
(4.58). 

Repeating the same argument leading to (4.58) one shows that on the event Mi, the event 

M2 := {^K+iit^+i) = C'K+,{t A+i) ^ -L + 2L1/3 , i = 1, . . . , K} , 

has probability at least {5{K + 1)^)^"*^. Thus, iterating, one concludes that 

P(^(l) = ^'(1)) ^ F{E) (4(K + l)2)-^+i {5{K + Ify^ ^ c(A)(5(i^ + 1))-^^ . 

This ends the proof of (4.55). D 

4.5. Proof of Theorem 1.2. Concerning the bound T^gj ^ L^"^^, we note that it follows from 
the same arguments used for the proof of Proposition 4.4. In fact, the situation is simpler here 
due to the absence of constraints on the number of crossings. We omit the details. 

Concerning the lower bound, we can actually prove that for some constant c > 0, one has 

T,^i^cL5/2(logL)-i. (4.61) 

We use an argument similar to that in (4.38). Recall that 

Tf., > ^ (4.02) 

for any / : 5 i— )■ M, where u is defined by (1.7) and 



nf,f)= X. ^[Var.J/)], 
xeOL 

Vx denoting the conditional probability v(^-\ay ^ y ^ x). Let C, denote the number of +1 in o", 
i.e. 

C(a)= Y, l(c7, = +1), C7G5. 

xeOL 
Note that C G {0, ...,-L}. Define the function f{a) = g{C{(j)/L), where g : [0,1] ^ [-1,1] 
is given by g{s) = —1 for s < 1/4, g{s) = 1 for s > 3/4, and by the linear interpolation for 
s G [1/4,3/4]. Since C = and Q = L have both positive probability uniformly in L, one has 
Var(/) ^ c for some constant c > 0. Let us estimate the Dirichlet form of /. We have 

V{t\ f) ^cL-^Yl '^(^) E P-(^)1(V4 ^ C ^ 3L/4) , 

for some constant c, where we write Px{(^) for the probability of a flip at x in a. Recall that the 
sum over x of the probabilities ^^(o") between consecutive crossings give a contribution of order 
1; this follows easily from Lemma 2.2. That yields '^xeo PxW) ^ ex for a suitable constant c, 
where x is the number of crossings (4.1). Therefore, adjusting the constant c: 

V{fJ)^cL-^v{x; L/4^C^3L/4). 

However, 

u{x, L/4^C^3L/A) ^u{x; X> clog L) + c (log L)iy(L/4^C^3L/4) . 



40 P. CAPUTO, H. LACOIN, F. MARTINELLI, F. SIMENHAUS, AND F.L. TONINELLI 



From Lemma 2.2 we deduce that i^ (x J X > c logL) = 0{L~'P) for some large p, provided c is 
large. On the other hand L/4 ^ C ^ 3L/4 implies that there exists a crossing in some position 
X G {—L + L/8, . . . ,L — L/8}. From the estimates of Lemma 2.2 this last event has probability 
0(L-i/2)^ so that 

u (L/4 ^ C ^ 3L/4) = 0{L-^/^) . (4.63) 

In conclusion, adjusting the constant c, one has !?(/, /) ^ cL^^'"^ logL. This ends the proof of 
(4.61). D 
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